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ABSTRACT 


A  statistical  analysis  of  the  output  signals  of  an  acousto-optic  spectrum  analyzer 
(AOS A)  is  performed  for  the  case  when  the  input  signal  is  a  continuous- wave  (CW). 

To  this  end,  a  statistical  model  of  these  output  signals  is  presented  as  a  basis  for  the 
different  analyses  along  with  some  numerical  algorithms  to  calculate  the  deterministic 
components  on  a  digital  computer. 

Using  this  model,  the  optimum  test  for  the  detection  of  a  known  frequency  is 
derived  and  its  performance  is  analyzed.  To  deal  with  the  unknown  frequency  situation, 
a  scheme  which  is  easily  implemented  with  a  finite  impulse  response  (FIR)  filter 
is  presented  and  its  performance  degradation  as  compared  to  the  case  of  a  known 
frequency  is  analyzed. 

The  frequency  estimation  problem  is  also  analyzed  and  the  Cramer-Rao  lower 
bound  on  the  variance  of  any  unbiased  estimator  is  calculated.  Since  the  Cramer-Rao 
bound  indicates  that  any  unbiased  estimator  would  exhibit  undesirable  characteristics, 
the  performance  of  the  peak-detecting  estimator  is  analyzed.  It  is  shown  that  this  latter 
estimator  is  biased  but  has  the  desired  characteristic  of  having  a  zero  average  bias. 

Finally,  the  power  estimation  problem  is  analyzed  under  the  assumption  that 
the  frequency  of  the  input  signal  is  known.  Under  this  assumption,  it  is  found  that 
the  maximum  likelihood  (ML)  estimator  is  an  efficient  estimator  which  means  that 
its  variance  is  the  lowest  possible  variance  of  any  unbiased  estimator.  The  effects  of 
inaccuracies  in  the  frequency  assumption  on  the  performance  of  the  ML  estimator  are 
also  analyzed. 


RESUME 


Une  analyse  statistique  des  signaux  generes  par  un  analyseur  de  spectre  acousto- 
optique  est  effectuee  pour  des  ondes  continues.  A  cette  fin,  on  presente  un  modele 
statistique  des  signaux  servant  de  base  aux  differentes  analyses  ainsi  que  des  algorithmes 
pour  le  calcul  des  signaux  en  l’absence  de  bruit. 

En  utilisant  ce  modele,  on  obtient  le  test  decisionnel  optimaie  pour  la  detection 
d’un  signal  dont  la  frequence  est  connue  et  on  analyse  sa  performance.  Pour  traiter 
du  cas  dont  la  frequence  est  inconnue,  on  presente  une  solution  qui  est  facilement 
executee  par  un  filtre  digital  dont  la  reponse  a  une  impulsion  est  finie  et  on  analyse  sa 
performance  degradee  par  rapport  au  cas  dont  la  frequence  est  connue. 

Le  probleme  de  l’estimation  de  la  frequence  est  egalement  analyse  et  la  limite 
Cramer-Rao  pour  la  variance  minimale  d’une  estimation  non-biaisee  est  calculee. 

Compte  tenu  des  characteristiques  indesirables  d’une  estimation  non-biaisee  mises  en 
lumiere  par  la  limite  Cramer-Rao,  on  calcule  la  performance  de  l’estimation  qui  consiste 
a  detecter  le  pic.  On  demontre  que  cette  derniere  estimation  est  biaisee  mais  qu’elle  a 
l’heureuse  characteristique  d’un  biais  moyen  nul. 


in 


En  dernier  lieu,  le  probleme  de  l’estimation  de  la  puissance  est  analyse  en  assumant 
que  la  frequence  du  signal  soit  connue.  Pour  ce  dernier  cas  ou  la  frequence  est  connue, 
on  demontre  que  l’estimation  a  probabilite  maximale  est  efficiente,  ce  qui  veut  dire  que 
la  variance  de  cette  estimation  non-biaisee  est  minimale  par  rapport  a  n’importe  laquelle 
des  estimations  non-biaisees.  Les  effets  de  1’inexactitude  de  notre  connaissance  de  la 
frequence  sur  la  performance  de  l’estimation  a  probabilite  maximale  sont  egalement 
analyses. 
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IV 


EXECUTIVE  SUMMARY 


DREO  is  investigating  the  use  of  Bragg  cell  receivers  for  the  development  of  Radar 
Electronic  Support  Measures  (RESM)  systems  for  the  armed  forces.  Perhaps  the  most 
mature  Bragg  cell  receiver  at  the  present  is  the  one-dimensional  configuration  which  has 
come  to  be  commonly  known  as  the  acousto-optic  spectrum  analyzer  (AOSA).  In  this 
report,  we  are  concerned  with  the  processing  of  the  output  signals  from  an  AOSA  for  the 
case  when  the  input  signal  is  a  continuous-wave  (C'W). 

To  this  end.  we  present  a  statistical  model  of  the  system  along  with  some  numerical 
algorithms  to  calculate  the  outputs  on  a  digital  computer.  Using  this  model,  we 
obtain  a  detection  algorithm  which  is  easily  implemented  and  we  also  characterize 
its  performance.  Next,  we  consider  the  frequency  estimation  problem  and  we  show 
that  an  estimator  of  any  significance  would  have  to  be  biased  especially  for  certain 
system  configurations.  In  order  to  provide  a  viable  solution  to  the  frequency  estimation 
problem,  we  propose  the  peak-detecting  algorithm  and  we  characterize  its  performance. 
Finally,  we  consider  the  power  estimation  problem  and  we  show  that  the  maximum- 
likelihood  (ML)  estimator  is  optimum  when  the  frequency  estimate  is  accurate  and  we 
also  analyze  the  effects  of  inaccuracies  in  the  frequency  estimate  on  the  performance  of 
•  the  ML  estimator. 

This  report  is  a  major  step  in  the  development  of  post-processing  algorithms  for 
'  the  AOSA.  Although  we  have  confined  our  analysis  to  the  case  of  CW  signals,  the 

derivations  could  be  applied  to  other  signal  types  and  system  configurations.  This 
report  contributes  to  a  better  understanding  of  the  critical  issues  and  tradeoffs  that 
are  involved  and  in  this  way  provides  a  basis  for  the  development  of  an  optimal  or 
suboptimal  solution  to  the  general  problem  which  is  more  complex. 
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1.0  INTRODUCTION 


1.1  Background 

Over  the  last  century,  electronic  systems  have  grown  to  predominance  in 
the  command,  control  and  communications  area  of  military  logistics.  Besides  the 
development  of  communications  and  radar  systems  for  military  applications,  the  desire 
to  disrupt  the  enemy’s  electronic  systems  has  fostered  the  development  of  a  new  area 
called  electronic  warfare  (EW).  The  purpose  of  EW  is  to  make  use  of  electromagnetic 
energy  to  determine,  exploit,  reduce  or  prevent  the  enemy’s  use  of  the  electromagnetic 
spectrum,  while  insuring  the  friendly  use  of  this  spectrum.  In  electronic  warfare 
applications,  receivers  can  be  used  to  intercept  signals  from  enemy  transmitters,  while 
jamming  transmitters  are  used  to  generate  false  information  or  noise  to  modify  the 
signal  received  by  the  enemy.  The  intercepted  information  could  be  used  to  identify 
and  evaluate  the  threat  associated  with  an  emitter  or  to  detect  the  parameters  of  the 
transmitter  which  would  help  in  a  jamming  operation.  In  this  context,  although  some 
information  on  the  transmitter  is  available,  the  transmitter  parameters  cannot  be  used 
in  the  design  of  the  receiver  since  the  radiating  sources  are  noncooperative.  Under 
such  conditions,  the  intercepting  receiver  must  be  designed  somewhat  independently 
of  the  transmitter  so  that  it  is  able  to  receive  a  variety  of  different  signals  and  extract 
the  desired  information.  For  that  reason,  the  design  of  such  receivers  is  complex  and 
requires  a  good  understanding  of  the  signal  environment  in  which  the  receiver  must 
operate. 

Microwave  electronic  warfare  receivers  for  this  kind  of  operation  can  be  divided 
into  the  following  groups  according  to  their  structures:  crystal  video  receivers, 
superheterodyne  (“superhet”)  receivers,  instantaneous  frequency  measurement  (IFM) 
receivers,  channelized  receivers,  compressive  (also  called  microscan)  receivers,  and  Bragg 
cell  receivers.  For  electronic  warfare  applications,  the  performance  of  the  latter  three 
types  of  receivers  is  expected  to  be  fax  better  than  that  of  the  former  three  types  in 
terms  of  the  width  of  the  spectrum  that  can  be  dealt  with,  the  signal  levels  that  could 
be  sucessfully  intercepted  and  the  number  of  signals  that  could  be  simultaneously  dealt 
with  by  the  receiver  [1,  pp.  3].  The  full  potential  of  these  receivers  has  not  yet  been 
reached  due  to  various  technical  difficulties  and  most  of  the  research  and  development 
presently  underway  in  this  area  is  aimed  at  overcoming  these  difficulties. 


1.2  Acousto-Optic  Spectrum  Analyzer 
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In  this  report,  we  will  solely  be  concerned  with  the  most  mature  Bragg  cell 
receiver  configuration  which  has  come  to  be  commonly  called  an  acousto-optic  spectrum 
analyzer  (AOSA).  In  this  receiver,  input  electrical  radio  frequency  (RF)  signals  axe 
first  transformed  into  spatial  patterns  which  modulate  a  light  beam,  generating  a 
spatial  distribution  of  the  light  intensity  which  is  sensed  by  a  photodetector  array.  This 
spatial  distribution  depends  on  the  parameters  of  the  signals  received  and  can  be  used 
to  discern  the  needed  information.  While  many  modulation  techniques  can  be  used 
in  optical  processors  (e.g.,  thermo-plastic  deformation  and  electro-optic  modulation), 
Bragg  cell  receivers  use  acousto-optic  modulation.  In  this  case,  the  electrical  signal  is 
converted  into  an  acoustic  wave  which  propagates  through  an  optically  transparent 
material  (Bragg  cell).  Through  the  elasto-optic  effect,  the  acoustic  wave  produces  a 
spatial  modulation  of  the  refractive  index  in  the  Bragg  cell.  When  a  coherent  light  wave 
is  passed  through  the  Bragg  cell,  the  refractive  index  modulation  (and  hence  the  electric 
signal  waveform)  is  impressed  onto  the  optical  wavefront  as  a  spatial  phase  modulation. 

A  suitable  optical  lens  system  is  used  to  convert  the  modulated  optical  wavefront  into 
a  spatial  intensity  modulation  corresponding  to  the  power  spectrum  of  the  Bragg  cell 
input  signal.  The  transformed  signal  is  read  and  converted  back  to  electrical  form  using 
a  linear  array  of  photodetectors.  Since  each  photodetector,  due  to  its  physical  size, 
provides  an  output  which  is  proportional  to  the  integral  of  the  spectrum  over  a  narrow 
portion,  we  can  describe  the  AOSA  as  a  form  of  channelized  receiver. 

Although  the  principles  discussed  above  have  been  known  for  many  years,  it  was 
only  fairly  recently  that  the  requisite  technologies  have  become  sufficiently  developed  to 
make  Bragg  cell  optical  processors  feasible.  Of  particular  importance  is  the  development 
of  lasers  with  adequate  output  levels  since  high  intensity  in  the  optical  wave  is  essential 
to  achieve  satisfactory  performance.  Other  significant  technological  advances  over  the 
last  decade  include  the  achievement  of  large  (>  100)  time-bandwidth  products  in  Bragg 
cells  and  the  development  of  large  photodetector  arrays.  While  much  work  remains  to  be 
done,  especially  in  the  areas  of  dynamic  range  and  output  rate  of  detector  arrays,  these 
recent  developments  in  optical  processor  technology  have  encouraged  the  exploration  of 
Bragg  cell  processors  for  microwave  receiver  applications. 

The  most  attractive  aspect  of  using  the  Bragg  cell  as  a  microwave  receiver  is  its 
potentially  extremely  small  size  and  low  cost.  Theoretically,  a  Bragg  cell  receiver  can 
perform  as  a  conventional  channelized  receiver  without  the  hundreds  of  filters  required  in 
such  a  receiver.  The  Bragg  cell  receiver  can  have  a  maximum  time-bandwidth  product 
of  approximately  1000  which  is  equivalent  to  a  channelized  receiver  with  1000  filters.  [1, 
pp.  150].  Furthermore,  the  development  of  integrated  optical  circuits  (IOCs)  makes  the 
integration  of  the  laser  source,  the  Bragg  cell  transducer,  output  detector  arrays,  and  the 
optical  lens  system  on  a  single  chip  possible.  An  integrated  optical  Bragg  cell  receiver 
could  have  a  volume  as  small  as  0.1  x  2  x  6  cm3  [2].  These  developments  make  the  Bragg 
cell  approach  the  most  attractive  electronic  warfare  receiver  for  airborne  applications  [1, 
pp.  150]. 
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1.3  Problem  Statement  and  Report  Organization 

In  this  report,  we  are  concerned  with  the  problem  of  processing  the  output  signals 
of  an  AOS  A  when  this  receiver  is  used  to  monitor  the  electromagnetic  environment  in 
an  effort  to  detect  the  presence  of  radar  signals  and  measure  their  respective  parameters. 
Such  E\\  receivers  must  be  able  to  cope  simultaneously  with  many  different  radar 
signals.  In  addition  to  radar,  other  types  of  signals  may  be  present  including  beacons 
and  transponders,  jammers,  missile  guidance  signals,  data  signals,  altimeter  signals, 
navigation  emissions,  and  identification  friend  or  foe  (IFF)  signals  [3,  pp.  1].  This  means 
that  there  is  a  large  variety  of  signals  that  we  can  expect  to  intercept.  Since  a  statistical 
analysis  dealing  with  ail  of  these  signals  at  once  would  be  extremely  complex  and  likely 
intractable,  we  will  limit  our  analysis  to  the  case  when  we  are  receiving  continuous- wave 
'  C'\\  )  signals.  In  the  case  of  a  CW  signal,  once  we  have  detected  its  presence,  there  are 
two  parameters  that  we  are  interested  in  estimating:  its  frequency  and  its  power.  We 
tei'l  that  this  is  a  good  place  to  start  as  it  will  contribute  to  a  better  understanding  of 
the  critical  issues  and  tradeoffs  and  provide  a  more  rigorous  and  systematic  foundation 
upon  which  optimal  or  suboptirnal  algorithms  can  be  derived  for  the  solution  of  the 
whole  problem.  We  will  therefore  address  the  problem  of  finding  efficient  algorithms  for 
the  processing  of  the  outputs  of  an  AOSA  in  order  to  fulfill  the  tasks  of  detecting  the 
presence  of  CW  signals  and  estimating  their  respective  frequency  and  power. 

To  begin  the  analysis  of  this  problem,  we  present  in  section  2.0  a  statistical  model 
of  the  output  signals  from  an  AOSA  which  will  form  the  basis  of  the  statistical  analyses 
performed  in  the  subsequent  sections.  Using  this  model,  we  then  consider  the  detection 
problem  in  section  3.0.  the  frequency  estimation  problem  in  section  4.0,  and  the  power 
estimation  problem  in  section  5.0.  Finally,  in  section  6.0  we  summarize  the  results  and 
discuss  the  needs  for  further  research. 


2.0  THE  SYSTEM  MODEL 


2.1  Introduction 

In  order  to  perform  a  'tat;  :  a  <d  analysis  of  the  performance  of  any  given  system, 
we  must  first  obtain  a  math*  m.ra  al  model  of  the  signals  and  the  noise  for  that  system. 
In  this  section,  we  pre-ent  a  m-Me]  for  the  signals  and  the  noise  at  the  output  of  an 
acousto-optic  spectrum  analy/.i  t  i  AOSA).  This  section  is  fundamental  to  this  report  as 
all  the  following  section^  will  be  based  on  it.  Indeed,  the  detection,  frequency  estimation 
and  power  estimation  analyse.-  performed  in  the  subsequent  sections  will  be  based  on  the 
model  presented  in  this  sect  a -n. 

In  the  next  subsection,  we  present  a  model  for  the  signals  at  the  output  of  an 
AOSA  assuming  the  noiseless  situation.  To  this  end,  we  present  the  AOSA  configuration 
that  will  be  analyzed  in  this  report  along  with  a  first-order  theory  of  operation.  Next, 
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we  present  a  mathematical  model  that  can  be  used  to  represent  the  signals  and  we 
give  some  numerical  algorithms  that  can  be  used  to  calculate  the  outputs  on  a  digital 
computer. 

In  the  third  subsection  of  this  section,  we  give  the  noise  model  that  will  be  used  in 
this  report  along  with  a  discussion  of  why  it  was  chosen.  Finally,  in  the  last  subsection 
of  this  section  we  conclude  by  giving  the  complete  signal  plus  noise  model  that  will  be 
used  in  the  subsequent  sections  for  the  different  statistical  analyses. 


2.2  Signal  Model 


2.2.1  AOSA  Configuration 

A  block  diagram  of  the  AOSA  configuration  of  interest  in  this  report  is  shown 
in  Figure  1.  The  first  component  of  this  system  is  the  laser  which  is  the  source  of  the 
optical  wave.  Since  the  beam  provided  by  the  laser  is  usually  relatively  small,  a  beam 
expander  is  required  so  that  a  plane  wave  commensurate  with  the  physical  size  of  the 
Bragg  cell  is  obtained.  This  is  usually  accomplished  with  a  series  of  lenses  which  are  also 
used  to  improve  the  quality  of  the  light  in  terms  of  getting  a  plane  optical  wave. 

The  key  component  of  this  system  is  the  Bragg  cell  which  acts  as  an  input  device 
by  transforming  the  input  electrical  signal  into  an  acoustic  wave  that  propagates  in  a 
transparent  medium  and  therefore  interacts  with  the  optical  wave.  The  phenomenon 
by  which  this  interaction  takes  place  is  called  acousto-optic  diffraction.  The  mechanism 
through  which  acousto-optic  diffraction  takes  place  is  due  to  the  fact  that  when  an 
acoustic  wave  propagates  in  a  transparent  medium,  it  induces  localized  refractive  index 
variations  via  the  elasto-optical  effect.  The  acoustic  wave  acts  like  a  moving  phase 
grating  which  may  diffract  portions  of  an  incident  light  beam  into  one  or  more  beams 
which  are  referred  to  as  the  different  diffracted  orders.  Provided  the  Bragg  cell  is  tilted 
by  the  proper  angle  6g  (the  Bragg  angle)  then  it  is  possible  to  obtain  onlv  one  diffracted 
order.  The  proper  value  for  9g  depends  on  the  frequency  of  the  optical  wave  and  the 
frequency  of  the  acoustic  wave.  (The  reader  in  search  of  an  intuitive  understanding 
of  acousto- optic  diffraction  can  consult  [4],  while  a  more  rigorous  derivation  from 
the  fundamental  Maxwell’s  equations  is  given  in  [5].)  An  important  characteristic  of 
acousto-optic  diffraction  is  that  the  diffracted  light  wave  will  be  phase  modulated  by  the 
phase  of  the  acoustic  wave  with  the  result  that  it  will  be  diffracted  at  an  angle  directly 
proportional  to  the  frequency  of  the  acoustic  wave.  The  Fourier  transform  lens  will 
then  map  the  twodimensional  Fourier  transform  of  this  diffracted  wave  onto  its  focal 
plane.  For  that  reason,  we  will  henceforth  refer  to  this  plane  as  the  frequency  plane  . 

(The  reader  who  desires  an  understanding  of  how  a  lens  can  perform  a  twodimensional 
Fourier  transform  can  consult  [6,  pp.  77-87].)  Since  we  assume  that  the  incident  optical 
wave  to  the  Bragg  cell  and  the  acoustic  wave  in  the  cell  are  both  plane  waves,  then  one 
dimension  of  the  frequency  plane  will  simply  be  the  one- dimensional  Fourier  transform 


I 
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of  the  input  signal  and  this  will  be  the  same  as  we  move  along  the  other  dimension  of 
the  frequency  plane. 

In  order  to  convert  the  result  of  this  optical  processing  to  electrical  form,  we 
place  a  linear  array  of  photodetectors  in  the  frequency  plane  (which  is  the  focal  plane 
of  the  Fourier  transform  lens).  These  photodetectors  perform  a  spatial  averaging  of  the 
light  intensity  which  represents  the  power  spectrum  of  the  input  signal.  In  this  report, 
we  will  also  assume  that  these  photodetectors  are  of  the  time-integrating  type.  This  is  a 
reasonable  assumption  as  many  practical  devices  operate  in  that  way,  which  means  that 
they  periodically  integrate,  sample  and  dump  the  light  intensity  (the  power  spectrum  of 
the  input  signal). 


2.2.2  Mathematical  Model 


From  the  standpoint  of  the  above  description  concerning  the  AOSA  configuration 
of  interest  in  this  report,  we  will  now  proceed  to  do  a  step  by  step  derivation  of  a 
mathematical  model  that  is  often  used  to  calculate  its  output  signals.  Using  this 
approach,  we  will  introduce  some  additional  practical  issues  and  show  how  these  are 
incorporated  in  the  model. 

As  a  starting  point,  if  u(t)  is  the  input  signal,  then  its  Fourier  transform  is 

roc 

w)  =  /  u(t)exp(—i2nft)dt. 


Now,  any  physical  Bragg  cell  would  have  a  finite  size  and  it  is  clear  that  the  light  in  the 
frequency  plane  can  only  represent  the  spectrum  analysis  of  that  portion  of  the  signal 
which  has  not  entirely  propagated  through  the  Bragg  cell  at  the  particular  point  in  time. 
In  effect,  this  means  that  the  AOSA  has  a  finite  time  aperture  over  which  to  perform 
spectrum  analysis  and  this  determines  its  fundamental  limit  in  spectral  resolution. 

To  account  for  this,  we  can  model  the  amplitude  and  phase  of  the  light  wave  in  the 
frequency  plane  as  a  sliding  window  spectrum  which  for  a  given  time  t  is 

/OO 

w(fl)u(t  —  /?)exp(— i27r//?)(f/9 

-OO 


where  w(0)  is  a  rectangular  function  to  account  for  the  fact  that  the  Bragg  cell  has  a 
finite  size  and  hence  the  input  signal  u(t)  will  appear  to  experience  truncation  in  time. 

It  turns  out  that  this  window  function  w{(i)  can  also  be  used  to  account  for 
some  other  practical  issues  concerning  the  Bragg  cell  and  the  laser  beam.  For  instance, 
the  window  function  can  be  used  to  account  for  the  fact  that  the  acoustic  wave  will  be 
attenuated  as  it  propagates  in  the  cell.  We  usually  account  for  this  by  changing  w(0) 
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from  a  rectangular  function  to  a  truncated  exponential  decaying  function.  The  extent 
to  which  it  will  be  important  to  account  for  this  attenuation  of  the  acoustic  wave  as  it 
propagates  through  the  cell  will  depend  on  the  Bragg  cell,  the  frequency  band  and  the 
specific  system  design  that  is  used. 


Still  another  practical  issue  that  the  window  function  u>(/?)  is  able  to  incorporate 
is  the  uniformity  of  illumination  on  the  Bragg  cell.  Most  practical  lasers  have  a 
Gaussian  shaped  profile  and  this  fact  can  be  used  to  obtain  a  windowed  spectrum 
analysis  with  lower  sidelobes.  This  effect  can  be  incorporated  in  the  model  through  the 
function  w(f3)  to  which  we  add  a  Gaussian  component.  The  extent  to  which  it  will  be 
important  to  account  for  this  practical  issue  will  depend  on  the  laser  and  the  specific 
system  design  that  is  used. 


Yet  another  component  of  our  AOSA  system  which  will  affect  our  mathematical 
model  is  the  linear  array  of  photodetectors  that  we  use.  First,  we  should  note  that  these 
photodetectors  are  sensitive  to  the  light  intensity  which  corresponds  to  the  squared 
magnitude  of  the  optical  wave.  Hence,  in  our  model  we  will  use  a  sliding  window  power 
spectrum  which  for  a  given  time  t  is 


/: 


w(/3)u(t  —  0)exp(—i2irfft)d0 


Secondly,  we  should  note  that  the  array  does  not  contain  an  infinite  number  of 
photodetectors  and  nor  are  these  of  infinitesimal  size.  Since  the  photodetectors  have 
a  certain  width,  they  will  spatially  integrate  the  light  intensity  which  corresponds  to 
a  frequency  integration  of  the  power  spectrum  of  the  input  signal.  We  can  therefore 
represent  the  output  of  an  individual  photodetector  as 


£ 


H(f-fk)  /  u>(/3)u(t  -  0)exp(-i2zf/3)d/3 


df 


where  H(f)  describes  the  spatial  response  of  that  photodetector.  Usually,  it  is  assumed 
that  H(f)  is  a  rectangular  function  and  the  bandwidth  that  it  represents  is  analogous 
to  the  bandwidth  of  a  filter  in  a  channelized  receiver.  However,  if  we  know  which  array 
of  photodetectors  will  be  used  and  there  is  some  data  available  concerning  the  profile  of 
t  ie  individual  detector  elements,  then  H(f)  can  be  used  to  take  this  into  account.  The 
frequency  that  corresponds  to  the  center  of  a  photodetector  element  is  called  and  we 
will  sometimes  refer  to  that  frequency  as  the  frequency  associated  with  the  fcth  detector. 


Thirdly,  as  we  mentioned  earlier,  we  will  assume  in  this  report  that  the 
photodetectors  are  of  the  time-integrating  type.  In  that  case,  we  get  that  the  ouput  of 
an  individual  photodetector  can  be  represented  as 


/•(>+!)/  1 

/  /  H(f-fk)  \  w(/3)u(t-  (3)exp(-i2nf0)dp\ 

J  jl  J—oo  \J-oo 


2 

df  dt 


-  a  - 


where  I  is  the  integration  time  of  the  photodetectors. 

To  summarize,  if  we  assume  that  an  integrating  photodetector  array  is  placed 
in  the  frequency  plane,  then  a  mathematical  model  that  can  be  used  to  describe  the 
signals  produced  by  the  AOS  A  and  the  one  that  will  be  used  throughout  this  report  is 
represented  by  the  following  equation 


r°°  r° 

Xjk=  /  /  HU-fk ) 

J  j  I  J  —  oo  J  —  c 


w(0)u(t  —  ,d)exp(-227r//3)d/? 


dfdt  (1) 


where  Xjj t  is  the  voltage  produced  at  the  output  of  the  Hh  detector  for  the  jth 
integrated  time  frame.  Basically,  this  equation  implies  that  the  instantaneous  light 
intensity  distribution  shining  on  the  array  of  photodetectors  is  the  magnitude  squared 
Fourier  transform  of  the  part  of  the  signals  that  are  contained  in  the  Bragg  cell  at  that 
time,  windowed  by  the  function  w(/3).  In  addition,  it  implies  that  each  photodetector  in 
the  frequency  plane  spatially  integrates  this  light  intensity  distribution  and  converts  it  to 
currents  which  are  integrated  and  sampled  at  periodic  intervals  of  duration  I. 

In  [7]  we  find  a  comparison  between  an  extended  version  of  this  model  and  some 
experimental  results,  which  serves  to  validate  our  model.  It  should  be  noted  that  the 
model  that  we  will  use  in  this  report  is  often  used  [8]  [9]  to  perform  analyses  on  the 
output  signals  of  an  AOSA. 

When  u(t)  is  a  pure  sinusoid  (which  corresponds  to  the  case  where  the  input 
signal  is  a  continuous-wave),  the  mathematical  model  represented  by  equation  (1)  can 
be  rewritten  in  a  more  convenient  form.  To  see  this,  we  note  that  equation  (1)  can  be 
rewritten  as 

Jz-O  +  l)/  roo  ra o  f  oo 

f  /  /  /  H(f-fk)w(a)w*(l3)u(t-a)u*(t-l3) 

]  I  J  —  oo  J  —oo  J  —oo 

exp\—i2irf(a  —  f3)}dad0dfdt.  (2) 


If 


u(t )  =  Acos(27r/ot  -f  <f>) 


(3) 


u(t  -  a )u*(t  -  /?)  =  A  cos  [2ir/0(t  -  a)  +  <f>]  ■  A  cos  [27r/0(t  -  (3)  +  <f>] 

=  ^{cos^tt/o*  -  27t/o(q  +  0)  4-  2<t>]  +  cos[27r/0(a  -  /?)]}. 


(4) 

(5) 


then 


Thus 


a2  f(j+l)I  roo  /■ oo  />oo 

X]k  =  -7T\  /  /  /  ^’(/“/*M«)u,*(^){cos[27r/o(a -^)] 

"  </  j/  J —oo  J  —oo  J  —  OO  k 

+  cos  [47t/0<  —  27r  fo(oi  +  fi)  +  2«ji]  |  exp  [— i2nf(ct  —  /?)]  (fa  d/?  d/ d£ 

.42 


(6) 


4  2  roo  roo  roo 

=  V  /  /  /  i/(/-/A:)«'(a)^*(^)exp[-227r/(o-^)] 

^  J  —oo  J  —oo  «/  — oo 

r(i+i)l  (  . 

j  |cos  [27r/o(a  —  /?)]  +  cos  [47r/0f  —  27r/0(o  +  /?)  +  2<^>j  j- da  d/? d/. 


(7) 


In  practical  situations  we  can  expect  that  /  >■  l/47r/0  (the  integration  time  of  the 
system  is  over  several  periods  of  the  CW  signal)  so  that 


exp  [— i27r/(a  -  0)\  J  jeos  [2i rf0(a  -  /3)\  +  cos  [47t/0<  -  27 r/0(a  +  ,3)  +  2<£]  j 

=  exp  [— i27r/(a  -  (3)\  {cos  [27r/0(a  -  /?)]  I 

f  sin  [4tt f0t  -  27r/o(or  +  ft)  +  2 <f>)  |0+1)/ 


dt 


4tt/0 


,  } 

V  -  I  J  1 

«  exp  [— i27r/(a  -  /?)]  cos  [27r/0(a  —  /?)]  J 

=  ^exp  [-i27t/(q  -  (3)]  |exp  [727t/0(q!  -  ft)]  +  exp  [-i27r/0(a  -  /?)]} 
=  ^{exP  H27r(Q  “  £)(/  ~  /o)]  +  exp  [-t'2ir(a  -  /?)(/  +  /o )] } - 

Substituting  (10)  in  (7)  we  get  that  for  /  1/47 r/0 

j  /-oo  COO  />oo 


(8) 

(9) 

(10) 


jk 


/OO  roo  roo  ( 

/  /  #(/  -  /fcMa)w*(/?){exp  [_i23r(a  -  /?)(/  -  /<,)] 

-oo  J  —  oo  J  —oo 

+  exp  [-*27t(q  -  /?)(/  +  f0)]}datdftdf.  (11) 


If  we  define 


/oo  |2 

w(t)  exp(— i27r/f)df 

-OO  I 

/oo  roo 

w(a)exp(—i2nfot)da  I  w*(ft)exp(i2irf  ft)dft 

-oo  J —oo 

/OO  roo 

/  w(a)w*(ft)  exp  [— i2wf(a  —  ft)]  do  d/? 

-oo  J  —oo 


(12) 


(13) 
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/OO  rOO 

/  w(a)w*(0)  exp  [-*27r(a  -  0)(f  +  /0)]  da < 

-oo  J  —  oo 


and  similarly 


/OO  /*00 

/  ^(o;)iy*(j9)  exp  [-i2ir(a  -  /?)(/  -  /0)]  da  dp. 

-oo  J  —  oo 


Substituting  (14)  and  (15)  in  (11)  we  get 


where 


A2 1 

X]k  &  T  /  #(/  —  /*)  [G(/  —  /o)  +  G(/  +  /o)]  d/. 

^  J  —  OO 

Now  if  we  assume  that  if(/)  is  symmetrical  about  /  =  0,  then 

f  HU  -  fk)G(f  +  fo)df  =  r  H(fk  ~  f)G(f  +  /o)d/ 

J  —  oo  J —  oo 

=  r  H  [/»-(/  +  /„)  +  /„]  G(f  +  f„)df 

J —oo 

=  r  BUk  +  h-  f'WW  =  nUk  +  /o), 

J  —oo 

H(f)  =  (°°  H(f  -  f')G(f')df 

J  —  oo 


which  is  the  convolution  between  the  functions  <?(/)  and  H(f).  Clearly  this  is  a 
function  centered  about  the  origin  with  bandwidth  equal  to  the  sum  of  the  bandwidths 
of  G(f)  and  H(f).  If  we  note  that 


«(/*  -  fo)  =  r  H(f  -  fk)G(f  -  fo)df 

J —oo 


then  by  substituting  (17)  and  (19)  in  (16)  we  get  that 


{H(fk  ~  fo)  +  H(fk  +  fo)}  ■ 


But  since  the  passband  of  a  practical  Bragg  cell  for  this  application  would  normally  be 
at  frequencies  which  are  high  with  respect  to  the  bandwidth  of  H(f),  then  H(fk  +  fo) 
will  be  very  small  compared  to  H{fk  —  fo),  and  thus  we  have  that  for  I  1/47 r/0 

A2 1 

Xit  a  —  HUk  -  ft). 


(21) 


(22) 
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We  shall  adopt  this  mathematical  model  throughout  this  report  and  use 

A2 1 

Xjk  =  — H(fk  -  /o), 


where  H(f)  is  the  convolution  between  the  functions  G(/)  and  H(f), 

H{f)  =  r  H(f  -  f')G(f')df 

J  —  OO 

with  G(f )  the  magnitude  squared  of  the  Fourier  transform  of  the  window  function 


(23) 


G(f)  =  f  w(t)exp(—i2-Kft)dt 

J  —OO 


(24) 


2.2.3  Numerical  Calculations 

When  the  input  signal  is  a  CW,  u(t)  —  Acos(2tt fQt  4-  <f> )  and  the  output  for 
the  fcth  photodetector  element  Xjk  can  be  obtained  by  appropriately  sampling  'H.(f) 
and  multiplying  by  a  scale  factor,  as  can  be  seen  from  equation  (22).  In  turn,  samples 
of  'H{f)  can  be  obtained  by  integrating  portions  of  G(f)  as  can  be  seen  from  equation 
(23).  If  H(f)  is  a  rectangular  function,  then  samples  of  H(f)  can  be  obtained  by  simply 
integrating  G(f)  over  the  appropriate  intervals  and  if  not,  then  we  need  to  window  G(f) 
according  to  H(f)  before  we  integrate. 

In  most  of  the  cases  of  interest,  there  is  no  closed  form  solution  to  equation 
(23)  and  hence  we  must  evaluate  'H(f)  numerically.  In  fact,  most  of  the  time  it  is  even 
difficult  to  obtain  a  closed  form  solution  for  G(f)  as  defined  in  equation  (24). 

The  derivations  performed  in  the  subsequent  sections  will  be  general  and  could 
be  applied  to  any  system  configuration.  However,  for  the  numerical  calculations  in  these 
sections  we  will  always  assume  that  H(f )  is  a  rectangular  function  of  unit  amplitude 
and  width  B  Hz,  symmetrical  about  /  =  0.  This  is  probably  a  reasonable  assumption  as 
it  is  unlikely  that  any  practical  H(f )  will  affect  our  calculations  to  any  great  extent.  In 
any  case,  the  calculations  could  easily  be  redone  for  any  H(f). 

In  addition,  the  numerical  calculations  performed  in  the  subsequent  sections  will 
only  consider  the  case  where  w(t)  is  a  rectangular  function  of  unit  amplitude  over  the 
interval  [0,r],  This  truncation  effect  due  to  the  finite  size  of  the  Bragg  cell  is  certainly 
the  most  important  factor  to  take  into  account.  However,  as  it  was  mentioned  in  the 
previous  subsection,  there  are  other  effects  which  may  or  may  not  need  to  be  taken 
into  account  depending  on  the  system  configuration  and  the  specific  components  that 
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are  used  in  the  system.  The  fact  is  that  there  could  be  many  possibilities  for  w(t), 
depending  on  the  application. 

Since  there  is  such  a  large  variation  in  t w(f),  in  this  report  we  have  chosen  to 
do  the  numerical  calculations  only  for  the  most  basic  case  where  w(t)  is  a  rectangular 
function.  We  assume  that  the  reader  interested  in  the  calculations  for  another  specific 
w(t)  will  redo  these  calculations  to  see  if  they  are  very  different  from  the  calculations 
performed  in  this  report.  However,  to  give  an  appreciation  of  the  effect  of  w(t)  on  the 
signals  produced  by  an  AOSA,  we  will  show  how  G(f)  can  be  calculated  for  some  of  the 
other  specific  <^ases  of  w(t)  and  how  this  will  affect  the  function  'H(f). 


2.2.3. 1  General  case 


A  family  of  window  functions  w(t)  that  is  found  to  be  useful  [10]  is  given  by 

,  .  2" 

w(t)  =  exp  j  -a(fQ)t  -  4 T2 


(H) 


,  *  1 
rect - - 

r  2 


(25) 


where  a  (in  nepers/sec)  accounts  for  the  acoustic  amplitude  attenuation  and  is  a 
function  of  the  input  signal  frequency,  T  accounts  for  the  profile  of  practical  laser  beams 
and  is  the  ratio  of  the  truncated  aperture  over  the  e~ 2  intensity  width  of  that  Gaussian 
profile,  and  r  (in  seconds)  is  the  truncated  aperture  which  is  related  to  the  physical 
length  of  the  Bragg  cell. 


If  a  =  0  (that  is,  if  we  assume  that  there  are  no  propagation  losses  in  the  Bragg 
cell),  we  have  that  u>(t)  is  a  truncated  Gaussian  profile.  To  see  the  effect  of  the  acoustic 
attenuation  on  that  window  function,  we  can  transform  equation  (25)  into  the  following 
form 


w(t)  =  exp 


t  1  or  \2  (ar)2 

7~2+8T*J  +  16T2 


(26) 


where  we  can  see  that  the  acoustic  attenuation  causes  a  shift  in  the  peak  position  of  the 
Gaussian  profile  as  well  as  a  decrease  in  the  peak  amplitude.  However,  the  general  shape 
of  the  window  function  is  preserved. 


In  general,  there  is  no  closed  form  solution  to  the  Fourier  transform  of  equation 
(25).  So  unless  some  simplification  or  approximation  is  done,  we  cannot  obtain  a  closed 
form  expression  for  G(f).  Once  we  have  obtained  G(f),  we  can  easily  obtain  'H(f) 
by  numerically  integrating  G(f)  over  finite  periods  since  we  assume  that  H(f )  is  a 
rectangular  weighting  function. 


We  can  estimate  G(f)  numerically  using  the  rectangular  rule  as 


G(f)  «  Ga(f)  = 


M -1 

At  ^  u;(mA<)  exp(— i'27r/mA<) 

m=0 


2 


(27) 
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If  we  let  At  =  t/M ,  where  M  is  the  number  of  points  from  w(t )  that  will  be  used  to 
calculate  one  sample  of  G(f),  then 


Ga(  f)  - 


A/-1 


(28) 


Defining  u  =  fr,  we  get 

Ga{v/r)  = 


or  if  we  express  the  complex  component  in  rectangular  coordinates,  we  get 

Ga(v/r )  =  |rtf(v)|2 


M-l 

T  V — ' 

/mr\  / 

'  —i2Trvm\ 

mD' 

m=0 

l  m)  exp  ( 

K  M  ) 

(29) 


where 


and  hence 


K(v) 


Em 


mJo1  h  (¥)  ™(2-r)  -  (¥) ™ 

M 

Ga(v/r )  =  r2  {[.Re(iir)]2  +  [/m(X)]2}  . 


(30) 


(31) 


(32) 


To  consider  the  error  associated  with  the  approximation  of  (32),  let 

M- 1 

Xa (/)  =  At  ^  w(mAt)exp(—i2nfmAt). 

m= 0 


(33) 


It  is  easy  to  see  from  (27)  that 


G.(f)  =  |X.(/)|J. 


(34) 


We  note  that  the  Fourier  transform  of  a  delayed  Dirac  delta  function  6(t  —  to)  is  e  ,2,r^fo. 


Thus  we  may  regard  (33)  as  the  Fourier  transform  of  At  Em=o  ^(mAt)6(t  —  mAt),  that 
is 

Xa(f )  =  T  |  At  ^  w(mAt)6(t  -  mAt)| .  (35) 


m= 0 


Since  u’(<)  is  zero  outside  the  interval  [0,r],  we  may  rewrite  (35)  as 

x,(f)  =  r{w(t)-At  £  e(t-mAt) 


(36) 


m  =  — oo 
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which  can  be  transformed,  using  the  well  known  property  for  the  Fourier  transform  of 
multiplied  time  functions,  to  the  following 


Xa(f)  =  W(f)  *  T  <  At  6(t~  rnAt) 


which  becomes 

xa(f)  =  w(f),  jr  <(/-£) 

m  =  — oo 


which  finally  becomes 

OO 

*«(/>  =  £  (37) 

m=  — oo 

So  we  see  that  the  approximation  of  equation  (32)  is  in  fact  an  aliased  version  of 
G(f).  We  also  see  that  Ga(f )  is  periodic  with  period  1/At,  hence  there  is  no  point  in 
calculating  Ga(f)  outside  the  range  of  frequencies 

-M  .  M 

~nZ~  <  f  <  7TZ 


since  At  =  r/M.  Or  if  we  use  the  normalized  version  Ga(v/r),  then  the  range  of 
frequencies  is 

-M  M 
___  <  v  <  — . 


In  summary,  we  have  shown  how  the  rectangular  rule  for  numerical  integration 
can  be  used  to  obtain  an  aliased  version  of  G(f).  The  resulting  normalized  numerical 
equation  is 


G.(v/t)  =  r3  {(fle(A:(»))]2  +  {Im(K(v))]2}  ,  (38) 


where  v  =  fr  and 


K(v) 


Sm=o  [*>  (ir)  ™  (£)  (2y )] 

M 


(39) 


where  M  is  the  number  of  sample  points  from  w(t)  that  will  be  used  to  calculate 
one  sample  of  G(f).  We  see  that  the  above  algorithm  is  closely  related  to  the  DFT 
algorithm  except  for  the  fact  that  it  can  be  used  to  obtain  samples  of  G(f)  at  any 
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frequency  instead  of  obtaining  those  for  a  fixed  set  of  frequencies.  Normalizing  w(t)  to 
the  following  equation 


u'(xt)  =  exp 


-L*-41«(z-  i) 


(40) 


where  L  =  or  and  .r  =  t/r,  makes  it  easy  to  calculate  equation  (38).  Ga(v/r)  is  periodic 
with  a  period  of  M  with  respect  to  e,  so  it  should  only  be  used  to  calculate  samples 
of  G(v/t)  for  the  range  <  v  <  4p.  Since  Ga(f)  is  an  aliased  version  of  G(/),  M 
should  be  made  large  enough  to  make  this  error  insignificant  and  we  should  not  attempt 
to  calculate  samples  of  G(f)  which  are  lower  than  a  certain  value.  The  value  of  M  and 
the  lowest  value  of  G(f)  that  we  attempt  to  calculate  using  Ga(v/r)  depend  on  G(/), 
so  determining  this  may  require  some  trial  and  error.  When  a  =  0  and  T  =  0,  w(t)  is 
simply  a  rectangular  window  and  G(f)  is  a  squared  Sine  function.  For  that  case,  the 
sidelobe  level  is  down  approximately  50  dB  from  the  peak  at  a  frequency  of  100/r.  For 
all  other  cases,  the  sidelobes  will  decrease  even  more  rapidly. 

Figure  2  shows  the  window  function  w(t)  for  the  case  of  ar  =  0.5,  T  —  1. 

Figure  3  shows  the  function  G(f)  calculated  from  equation  (38)  using  M  =  200  for  the 
same  case  of  ar  =  0.5.  T  —  1.  Finally,  Figure  4  shows  Ff(/)  for  that  same  case  assuming 
the  photodetectors  have  a  bandwidth  of  B  =  1/r.  It  should  be  noted  that  7d(/)  is 
symmetrical  about  /  =  0  even  though  this  is  not  shown  on  Figure  4.  This  last  figure  is 
the  result  of  numerically  integrating  the  function  G(f).  For  that  reason  it  is  useful  to 
have  a  numerical  algorithm  that  can  calculate  samples  of  G(  f)  at  any  frequency  because 
most  numerical  integration  programs  require  a  function  which  can  do  that. 

2. 2. 3. 2  Untruncated  Gaussian  windowing 

We  have  showed  how  the  output  of  the  AOSA  could  be  calculated  using  an 
algorithm  which  has  calculations  similar  to  those  encountered  in  the  direct  evaluation  of 
the  DFT.  But  it  has  long  been  recognized  that  the  DFT  is  computationally  expensive. 

In  fact  it  is  presumed  that  Carl  Friedrich  Gauss,  the  eminent  German  mathematician, 
developed  an  algorithm  that  could  simplify  its  calculation  as  early  as  the  year  1805  [11]. 

It  is  therefore  worthwhile  to  note  that  if  we  assume  a  =  0  and  the  Gaussian 
shaping  is  not  truncated  by  the  physical  size  of  the  Bragg  cell,  then  we  can  get  a  closed 
form  expression  for  G(f).  Harms  and  Hummels  have  used  this  result  in  [9]  and  have 
calculated  the  transform  of  w(t)  by  means  of  a  line  integral  in  the  complex  plane. 

In  this  case  we  have  that 


it'(t)  =  exp 


-4T2 


(41) 
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FIGURE  4:  H(f)  FOR  ar  =  0.5,  T  =  1,  tB  =  1 


so  that 

.  2  7T  /  — 7T2U2\ 

G(t-/r)  =  r  —  exp  (-jyj- J  (42) 


where  again  u  =  /r.  This  greatly  simplifies  the  calculation  of  H(/)  but  it  should  only  be 
used  when  T  is  large.  (In  [9]  this  approximation  for  T  =  1.63  is  used  and  it  is  claimed 
that  the  resultant  error  is  negligible.) 

Figure  5  shows  the  window  function,  w(t)  for  the  case  of  T  =  2  while  Figure  6 
shows  the  function  G(f  )  calculated  using  equation  (42).  Finally,  Figure  7  shows  7i(f) 
for  that  same  case  assuming  the  photodetectors  have  a  bandwidth  of  B  =  1/r. 

2. 2. 3. 3  Truncated  exponential  windowing 

If  we  assume  that  T  =  0,  that  is  if  we  assume  that  the  laser  has  no  Gaussian 
shaping  but  we  still  want  to  take  into  account  the  fact  that  there  is  attenuation  of  the 
acoustic  wave  as  it  propagates  through  the  Bragg  cell,  then  w(t)  is  given  by 


w(t)  =  exp(— af)rect 


(43) 


dB  RELATIVE  TO 


8 


-  19  - 


FIGURE  7:  H{f)  FOR  ar  =  0,  T  =  2,  tB  =  1 


and  so 


G(v/t) 


2  [l  -  2e  L  cos(27rv)  +  e  2L] 
L2  +  (27ru)2 


(44) 


Figure  8  shows  the  window  function  w(t)  for  the  case  of  ar  =  0.5  while  Figure  9 
shows  the  function  G(f )  calculated  using  equation  (44).  Finally,  Figure  10  shows  H(f) 
for  that  same  case  assuming  the  photodetectors  have  a  bandwidth  of  B  =  1/r. 


2. 2. 3. 4  Rectangular  windowing 

As  it  was  mentioned  earlier,  if  we  ignore  the  acoustic  loss  of  the  Bragg  cell  and 
the  Gaussian  shaping  of  the  laser,  then  w(t)  is  simply  a  rectangular  window  whose 
duration  is  determined  by  the  physical  size  of  the  crystal  or  the  width  of  the  light  wave 
impinging  on  the  cell.  For  that  case  it  is  easy  to  show  that  w(t)  is  as  shown  in  Figure  11 
and  G(f )  is  given  by  the  following  equation 


G(f)  =  t2 


sin2(27r/r/2) 

(2nfr/2)2 


(45) 


which  is  shown  in  Figure  12.  The  resulting  function  H(f)  for  tB  =  1  is  shown  in 
Figure  13. 
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FIGURE  8:  WINDOW  w(t)  WHEN  ar  =  0.5,  T  =  0 


FIGURE  9:  G(f)  WHEN  ar  =  0.5,  T  =  0 


LINEAR  AMPLITUDE  dB  RELATIVE  TO  i 


FIGURE  10:  H{f)  FOR  or  =  0.5,  T  =  0,  tB  =  1 


TIME 


FIGURE  11:  WINDOW  w(t)  WHEN  or  =  0,  T  =  0 
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2.3  Noise  Model 

In  the  previous  subsection,  we  presented  the  AOSA  configuration  of  interest  in  this 
report  along  with  a  mathematical  model  to  calculate  the  deterministic  component  of 
the  photodetector  output.  In  addition  to  the  deterministic  component  due  to  the  input 
signal,  the  photodetector  output  is  corrupted  by  noise  from  a  variety  of  sources.  Among 
the  possible  noise  contributors  we  can  include  the  laser,  the  microwave  front-end,  the 
photodetectors,  vibration  and  scattered  light. 

In  the  application  that  we  are  considering,  we  are  always  striving  to  design 
receivers  with  high  dynamic  range  and  good  sensitivity.  When  designing  an  AOSA  for 
this  application,  we  usually  find  that  the  photodetectors  are  the  bottleneck  in  terms  of 
those  requirements.  This  means  that  the  photodetector  noise  is  dominant  among  the 
other  sources  of  noise. 

In  [12]  we  performed  experimental  measurements  on  the  noise  of  an  avalanche 
photodiode  (APD)  array.  This  array  is  of  special  interest  for  the  design  of  systems 
to  monitor  the  electromagnetic  environment  because  of  its  high  sensitivity.  The 
measurements  were  done  using  an  analog  amplifier  which  would  typically  be  used  for 
this  application.  It  was  found  that  the  distributions  governing  the  noise  present  in  the 
photodetector  element  outputs  were  Gaussian  or  normal  distributions.  The  means  of  the 
probability  density  functions  were  found  to  vary  significantly  from  element  to  element 
although  their  variances  were  roughly  the  same. 

In  light  of  these  findings,  we  will  assume  in  this  report  that  the  signal  components 
of  the  detector  elements  are  corrupted  by  additive  Gaussian  random  variables 
independent  of  each  other  but  having  identical  variances.  In  addition,  we  will  assume 
that  these  random  variables  have  zero  means  since  in  practice  these  noise  offsets  would 
most  likely  be  subtracted  out.  Hence,  the  noise  model  that  will  be  used  in  this  report  is 
that  each  photodetector  output  will  be  corrupted  by  independent,  identically-distributed 
zero-mean  Gaussian  variables  with  variance  a2. 


2.4  Signal  Plus  Noise  Model 

Having  presented  the  signal  model  and  the  noise  model  that  will  be  used  in  this 
report,  we  will  now  present  the  complete  signal  plus  noise  model  which  will  be  used  in 
the  subsequent  statistical  analyses.  To  this  end,  let  us  first  define 

R  =  {ri>r2>r3, . . .  ,r/v}  (46) 


as  the  received  vector  which  is  used  to  represent  the  photodetector  outputs  after  any 
integration  time  frame.  We  will  sometimes  refer  to  the  photodetector  outputs  as  the 
pixel  outputs  since  this  is  an  expression  which  is  often  used  in  this  field. 
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We  note  that  R  is  a  vector  in  an  jV-dimensional  space  where  N  is  the  number  of 
photodetector  elements  used  in  the  linear  array.  We  also  note  that 

rt=mi(f0,A)  +  ni,  i  =  1,2,3, ...  ,JV  (47) 

were  the  n,’s  are  zero-mean,  independent  identically-distributed  random  variables 
with  variance  a2.  It  is  important  to  realize  that  nif.  =  Xjk  in  equation  (1).  We  have 
chosen  to  henceforth  use  the  later  notation  (i.e.  to  reflect  the  fact  that  the  signal 
components  of  the  r,  ’s  are  actually  the  means  of  these  random  variables.  It  is  equally 
important  to  note  that  the  m/s  are  completely  defined  once  we  know  the  frequency  fQ 
and  the  amplitude  A  of  the  input  signal. 

In  the  following  sections,  we  will  use  this  signal  plus  noise  model  and  proceed  to 
perform  statistical  analyses  in  an  effort  to  obtain  efficient  post-processing  algorithms. 

We  will  look  at  the  detection  as  well  as  the  frequency  and  power  estimation  problems. 

3.0  DETECTION 

3.1  Introduction 

The  first  problem  that  we  face  when  processing  the  output  signals  of  an  AOSA  is 
the  detection  problem.  That  is,  we  must  first  make  a  decision  as  to  whether  one  or  more 
signals  are  present  before  we  try  to  estimate  their  respective  frequency  and  power. 

In  order  to  find  an  optimal  solution  for  this  problem,  we  must  first  define  a  criteria 
upon  which  this  optimality  is  to  be  decided.  One  criterion  that  is  commonly  used  in 
classical  detection  theory  is  the  Bayes’  criterion.  The  use  of  this  criterion  requires  the 
existence  of  a  priori  statistics  of  the  observed  signal.  This  criterion  cannot  be  used  in 
the  present  circumstance  since  no  meaningful  a  priori  statistics  can  be  assigned.  The 
standard  method  of  dealing  with  such  situations  is  to  use  the  Neyman-Pearson  criterion. 

The  Neyman-Pearson  criterion  considers  the  conditional  probability  of  deciding 
that  a  signal  is  present  given  that  there  is  in  reality  no  signal  present  (the  false  alarm 
■probability  Pp )  and  the  conditional  probability  of  deciding  that  a  signal  is  present  given 
that  there  is  indeed  a  signal  present  (the  probability  of  detection  Pd).  In  applying  the 
Neyman-Pearson  criteria,  we  like  to  design  a  test  which  minimizes  Pp  and  maintains  Pd 
as  large  as  possible.  These  turn  out  to  be  two  conflicting  objectives  and  so  a  tradeoff 
must  be  made.  A  specific  aim  in  applying  the  Neyman-Pearson  criterion  is  to  maximize 
Pd  while  maintaining  Pp  less  than  a  specified  amount  £,  that  is  maximize  Pd  subject  to 
Pf  <  £.  This  is  the  criterion  that  we  will  use  in  the  remaining  subsections  of  this  section 
in  dealing  with  the  detection  problem  we  face. 

The  first  thing  to  note  about  this  detection  problem  is  that  we  are  dealing  with 
a  composite  hypothesis  situation.  This  is  because  we  seek  to  detect  the  presence  of  a 
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signal  of  which  we  do  not  know  the  frequency  or  the  amplitude.  We  sometimes  refer  to 
these  unknown  parameters  of  a  composite  hypothesis  problem  as  unwanted  parameters 
since  we  do  not  care  about  them  for  the  purposes  of  detection;  we  simply  want  to  decide 
whether  one  or  more  signals  are  present  irrespective  of  these  parameters. 

In  the  next  subsection,  we  simply  assume  that  these  parameters  (the  frequency 
and  the  amplitude)  are  known  and  we  go  about  deriving  the  optimum  test  and  its 
performance.  This  is  a  logical  thing  to  do  because  in  certain  instances,  by  assuming  the 
unwanted  parameters  are  known,  we  obtain  a  test  which  turns  out  to  be  independent 
of  these  parameters.  We  refer  to  such  a  test  as  a  uniformly  most  powerful  (UMP)  test. 

It  turns  out.  as  we  will  show  in  the  next  subsection,  that  the  optimum  test  we  obtain 
is  in  fact  independent  of  the  amplitude  of  the  input  signal.  However,  it  is  dependent 
on  its  frequency.  In  the  third  subsection,  we  deal  with  this  by  using  the  test  obtained 
in  section  3.2  to  decide  on  the  presence  of  a  number  of  possible  frequencies  within  the 
bandwidth  of  our  receiver.  We  obtain  a  scheme  which  has  a  good  performance  and 
which  is  easily  implemented. 

3.2  Detection  of  a  Known  Frequency 

In  this  subsection  we  will  find  the  optimum  detection  test  under  the  Neyman- 
Pearson  criterion  for  the  situation  where  we  want  to  decide  whether  a  CW  signal  with  a 
given  frequency  is  present  or  not.  To  this  end,  we  let 

R  =  {rl  1  r2>  r3>  •  •  •  >  rN  } 

be  the  received  vector  where  the  r,’s  are  the  pixel  output  values  for  a  given  frame.  We 
can  consider  this  problem  as  a  binary  hypothesis  test  where  Hq  is  the  hypothesis  that  no 
signal  is  present  and  Hi  is  the  hypothesis  that  a  signal  is  present  and  hence  we  can  write 

H0  :  r,  =  n„ 

Hi  :  r,  =  m,(/0,.4)  +  n,,  i  =  1,2,3, ...  ,N 

where  the  n,’s  are  zero-mean,  independent  identically-distributed  Gaussian  random 
variables  with  variance  cr2,  the  m,'s  axe  the  signal  components,  N  is  the  number  of 
pixels  in  the  photodetector  array  and  /0,  A  are  respectively  the  frequency  and  amplitude 
of  the  input  signal  as  defined  in  equation  (3).  Referring  to  equation  (1),  we  have  that 
772*  =  Xjk,  but  we  will  henceforth  use  this  latter  notation  to  denote  the  fact  that  the 
signal  components  are  the  means  of  the  Gaussian  random  variables  at  the  output  of  the 
AOS  A.  From  (22)  we  have  that 

4  2/ 

m,  =  — —H,  (49) 

4 

when-  it  is  understood  that  the  7d,’s  are  the  samples  of  ?f(/)  corresponding  to  the 
frequency  and  pixel  under  consideration  (i.e.  H,  =  'H{ft  -  /o))- 


It  is  well  known  [13,  pp.  34]  that  the  optimum  test  according  to  the  Neyman- 
Pearson  criterion  is  satisfied  by  the  likelihood  ratio  test 


A  (R)  >  X 


where  the  likelihood  ratio  test  is  defined  as 


A(R)  -  mw) 


where  p(R  |  Ht)  is  the  conditional  probability  density  function  of  the  received  vector  R 
given  that  hypothesis  Ht  is  true  and  A  is  some  threshold  depending  on  the  constraint 
Pf  <  £•  This  latter  quantity  can  be  determined  from  noting  that 


PF  =  J™  p(A(R)\H0)dA(R). 


It  is  easy  to  show  that  for  our  problem  as  it  is  defined  in  (48),  the  above 
conditional  probabilities  are  as  follows 


,5  ,  rr  x  TT  1  [~(r»  “  m«)2 

-g>-  j 


and  hence  the  likelihood  function  can  be  written  as 


A  (R)  =  exp 


=  exp 


2 


But  because  the  natural  logarithm  is  a  monotonic  function  and  both  sides  of 
equation  (50)  must  be  positive,  the  likehood  ratio  test  is  equivalent  to  the  test 


ln{A  (R)} 


ln(A) 


(57) 
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which  is  very  useful  given  the  form  of  A (R)  as  shown  in  equation  (56). 


Indeed,  we  can  say  that  a  test  satisfying  the  Ney man- Pearson  criterion  for  our 
problem  is  as  follows 


2<r2 


H  i 

J  ln(A). 
Ho 


(58) 


This  is  equivalent  to  the  test 


N 

E 

i=i 


r,m, 


H  i 
> 
< 
Ho 


7, 


(59) 


where  -y  is  some  constant. 

It  should  be  noted  that  hypothesis  H i  is  a  composite  hypothesis  because  it  contains 
the  unwanted  parameter  A  and  the  test  of  equation  (59)  cannot  be  used  unless  we  know 
the  value  of  A.  However,  since 

A2l-u 

mi  =  ——Hi 


we  can  reduce  the  above  test  to  the  following 


H  Hi 

i=l  Ho 


(60) 


because  A2 1  is  a  positive  quantity.  If  we  let 

N 

Y  =  J2r'n' 

1=1 


then  the  above  test  becomes 


Hi 

> 

< 

Ho 


7 


and 

Pf=  rP(Y\H0)dy. 

Jy' 


(61) 


(62) 


The  above  test  is  what  is  called  a  uniformly  most  powerful  (UMP)  test  because 
it  is  completely  defined  (including  threshold)  for  a  given  Pp  without  knowledge  of  the 
parameter  A.  which  is  the  amplitude  of  the  input  signal.  Of  course,  the  performance  of 
this  test  will  l>e  a  function  of  this  parameter  as  will  be  shown  shortly. 
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In  order  to  find  the  performance  of  this  test,  we  note  that  if  H0  is  true  then  Y  is  a 
Gaussian  random  variable  with  zero  mean  and  variance 


N 


And  similarly,  if  H\  is  true  then  Y  is  a  Gaussian  random  variable  with  mean 

N 

rriil-ii.  =  - 

4 


,2  r  N 


Y  =  J2rntHt  =  ±±J2'H? 

1=1  t=i 


and  variance 


N 


From  this  we  have  that 


Pf 


dx 


or 


where 


erfc(t)  =  ~^= 
V7r 


(63) 


(64) 


(65) 


and 


,D  r_i 

J-yl  \JlWOY 


exp 


2<7y 


dx 


(66) 


or 


(67) 


Figure  14  shows  the  performance  of  the  test  of  equation  (60)  as  a  function  of  the 
parameter 

N 


FIGURE  14:  PERFORMANCE  OF  THE  MATCHED  FILTER 
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which  serves  as  a  figure  of  merit.  Indeed  z  is  the  S/N  ratio  times  the  integration  time  of 
the  photodetectors  (I)  times  a  factor  which  quantifies  the  performance  gain  due  to  the 
matched  filter.  We  can  see  that  the  performance  gain  due  to  the  matched  filter  depends 
on  the  number  of  pixels  that  are  used  in  the  algorithm.  In  Appendix  I  we  calculate 
=  i  for  different  values  of  the  constant  tD.  Of  course,  in  practice  we  would  not 
use  .V  values  of  the  7t,'s  in  equation  (GO)  to  implement  this  test  because  some  of  these 
values  will  be  insignificantly  small.  Using  Figure  14  and  equation  (68)  we  could  find  a 
reasonable  number  of  Ti,' s  to  use  in  our  test  for  a  given  system  configuration. 

3.3  Detection  of  an  Unknown  Frequency 

In  the  previous  subsection,  we  looked  at  the  optimum  detection  test  using  the 
Xeyman-Pearson  criterion  for  the  case  of  a  signal  with  known  frequency.  We  have  found 
that  the  optimum  test  for  that  case  is  the  matched  filter  and  we  have  also  characterized 
its  performance.  However,  for  our  application  which  is  to  use  an  AOSA  to  monitor  the 
electromagnetic  environment,  we  are  not  only  interested  in  testing  for  the  presence  of 
one  known  frequency,  but  rather  we  are  interested  in  testing  for  the  presence  of  a  range 
of  unknown  frequencies. 

We  can  approach  this  problem  by  dividing  the  bandwidth  of  our  receiver  into  a  set 
of  possible  frequencies  such  that  the  resultant  sequence  forms  a  fine  mesh  over  the  whole 
bandwidth.  Having  done  this,  we  can  then  use  the  matched  filter  to  sequentially  test  for 
the  presence  of  each  one  of  these  possible  frequencies.  We  know  that  this  test  will  be 
optimum  if  the  input  signal  frequency  corresponds  to  one  of  the  possible  frequencies  that 
we  have  selected  for  our  set. 

A  natural  and  convenient  spacing  for  this  set  of  possible  frequencies  is  B,  the 
frequency  width  of  each  photodetector.  The  reason  for  this  is  that,  once  we  have 
selected  the  length  of  the  matched  filter,  let  us  say  n  pixels,  then  we  can  simply  slide 
the  matched  filter  over  the  whole  array,  comparing  its  output  against  a  threshold.  This 
is  simply  a  Finite  Impulse  Response  (FIR)  filter  algorithm  and  is  easily  implemented  as 
photodetector  arrays  often  use  serial  output  structures. 

The  advantage  of  this  scheme  is  that  we  need  not  change  the  threshold  value  that 
we  have  selected  according  to  the  false  alarm  probability  that  we  are  willing  to  live  with. 
Of  course,  something  different  will  have  to  be  done  for  the  frequencies  near  the  edge  of 
the  array.  We  could  simply  not  consider  the  edge  frequencies  if  the  value  of  n  is  small 
or  we  could  change  the  value  of  the  threshold  for  these  edge  frequencies.  When  a  strong 
signal  is  present,  the  output  of  the  FIR  filter  will  likely  exceed  the  threshold  for  more 
than  one  frequency,  but  this  is  acceptable  since  we  are  using  this  scheme  to  detect  the 
presence  of  signals  and  not  to  estimate  the  frequencies  of  those  signals. 

This  scheme  will  optimally  detect  the  presence  of  all  the  frequencies  selected  in 
our  set  of  possible  frequencies  under  the  constraint  of  a  given  false  alarm  probability. 

Wo  might  want  to  know  how  this  detection  scheme  performs  for  the  frequencies  that 


are  not  included  in  this  set;  we  cannot,  after  all  expect  to  have  a  frequency  line  up 
with  a  photodetector  center  frequency.  To  answer  this  question,  let  us  consider  the 
more  general  case  of  how  the  performance  of  the  test  of  equation  (60)  is  affected  when 
the  used  in  the  test  axe  different  from  the  Tii's  that  correspond  to  the  signal. 

This  could  happen  when  the  frequency  of  the  input  signal  is  different  from  the  one  we 
assumed  with  the  Ti,'s  of  the  filter  or  when  the  window  function  w(t)  that  we  used  to 
calculate  the  Ti,'s  was  not  accurately  measured  or  for  any  other  reason. 

To  calculate  this  effect  on  performance  let 


Y,r*nx 


H , 
Ho 


be  the  test,  where  the  Tit's  are  used  in  the  test  but  do  not  actually  correspond  to  the 
Tit's  of  the  signal.  In  this  case,  if  Hq  is  true,  Y  is  a  Gaussian  random  variable  with  zero 
mean  and  variance 


a2?  =  a2Y,Hl 


But  if  Hi  is  true,  then  is  a  random  variable  with  mean 

N  A2  T  N 

E{V-}  =£*>.-  =  — 


and  variance 


WE*?- 


It  is  easy  to  show  that  the  performance  of  the  test  of  equation  (69)  is  the  same  as  the 
performance  of  the  UMP  test  as  shown  in  Figure  14  except  that  for  this  case 

<’v/E,=. «?  ' 


We  can  therefore  take  the  ratio  of  equation  (70)  anr*  (68)  to  find  the  degradation 
in  defection  performance  when  the  input  signal  frequency  is  between  two  pixels  but 
the  matched  filter  used  is  designed  to  test  for  the  case  when  the  input  signal  frequency 
corresponds  to  one  pixel.  Figure  15  shows  this  degradation  for  different  values  of  the 
constant.  tD  assuming  rectangular  windowing.  We  have  used  a  matched  filter  with 
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25  taps  to  calculate  the  curves  of  Figure  15  and  this  assures  us  that  all  the  7{,'s  that 
were  dropped  were  insignificantly  small.  In  Figure  15,  we  have  only  calculated  the 
degradation  over  one  period  since  the  degradation  is  periodic  with  period  B.  As  an 
illustration  of  how  to  interpret  Figure  15,  a  degradation  of  1  dB  means  that  the  input 
signal  power  would  have  to  be  1  dB  higher  to  obtain  the  same  probability  of  detection 
as  if  the  input  signal  frequency  corresponded  to  the  frequency  of  a  pixel,  everything  else 
being  the  same. 

We  can  see  from  this  figure  that  for  a  given  tB,  the  worst  case  degradation  of  the 
matched  filter  occurs  when  the  input  signal  frequency  is  exactly  in  the  middle  of  two 
fk's.  We  also  see  that  this  worst  case  degradation  increases  for  larger  values  of  tB. 
However,  we  should  note  that  even  for  tB  =  3  the  worst  case  degradation  is  less  than 
1.5  dB,  a  relatively  small  figure.  Another  important  fact  to  note  about  Figure  15  is  that 
most  of  the  degradation  is  concentrated  in  a  relatively  small  frequency  range.  This  can 
be  seen  by  the  steepness  of  the  curves  at  the  point  where  the  input  signal  frequency  is 
exactly  in  the  middle  of  two  fk's. 

Hence  we  can  see  that  the  sliding  matched  filter  scheme  has  a  very  good 
performance  for  the  detection  of  unknown  frequencies.  The  worst  case  degradation  is 
relatively  small  and  most  of  the  degradation  occurs  over  a  relatively  small  frequency 
range. 
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4.0  ESTIMATION  OF  THE  FREQUENCY 


4.1  Introduction 

The  most  important  function  of  the  AOSA  receiver  is  to  determine  the  frequency  of 
the  input  signals.  In  this  section,  we  will  consider  the  problem  of  processing  the  output 
signals  of  the  AOSA  for  the  purpose  of  estimating  the  input  signal  frequency.  We  will 
approach  this  problem  using  classical  estimation  theory  as  expounded  by  Van  Trees  [13]. 

Classical  estimation  theory  is  divided  into  two  main  branches:  the  first  is  random 
parameters  estimation  and  the  second  is  nonrandom  parameters  estimation.  Random 
parameters  estimation  requires  some  knowledge  of  the  a  priori  probabilities  for  the 
outcome  of  the  values  of  the  parameters  that  we  wish  to  estimate.  In  many  cases, 
and  certainly  in  the  case  of  a  receiving  system  that  monitors  the  electromagnetic 
environment,  it  is  unrealistic  to  treat  the  unknown  parameter  as  a  random  variable. 
When  dealing  with  a  nonrandom  parameter  situation,  we  must  work  with  the  a 
posteriori  density  functions  and  evaluate  the  performance  of  different  estimators  by 
considering  their  bias  and  variance. 

A  useful  landmark  for  the  nonrandom  parameter  situation  is  the  Cramer-Rao 
bound  which  gives  us  the  lower  bound  on  the  variance  of  any  unbiased  estimator.  In 
the  next  subsection,  we  present  this  bound  and  show  its  dependence  on  basic  system 
design  parameters.  We  shall  see  that  this  bound  is  useful  since  it  shows  us  that,  for  this 
problem,  any  practical  estimator  of  any  significance  would  have  to  be  biased.  Hence  in 
the  third  subsection  of  this  section  we  present  the  performance  of  a  practical  estimator, 
the  peak-detecting  estimator.  This  estimator  is  a  biased  estimator  although  it  has  the 
desired  characteristic  of  having  a  zero  average  bias. 

The  derivations  in  this  section  will  be  kept  as  general  as  possible  and  could  be 
applied  to  any  system  independent  of  the  specific  window  function  used.  However, 
the  numerical  calculations  will  be  performed  for  the  case  of  the  rectangular  window 
function.  These  calculations  could  easily  be  done  for  any  window  function  if  we  were 
given  the  specifications  of  a  particular  system,  but  we  have  chosen  to  do  them  only 
for  the  rectangular  case  to  minimize  the  number  of  computations  which  are  already 
relatively  extensive  even  for  this  latter  case. 
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4.2  Cramer- Rao  Bound 

Using  the  model  described  in  section  2.0,  we  will  in  this  subsection  present  the 
Cramer-Rao  bound  for  the  estimation  of  the  input  signal  frequency.  To  this  end,  we  let 

R  =  { r l '  r 2 ,  r - rN} 


be  the  received  vector  where  the  r,'s  are  the  pixel  output  values  for  a  given  integrated 
time  frame.  We  can  write 

r,  =  +  nit  1  =  1,2,3 - ..V  (71) 

where  the  n,'s  are  zero-mean,  independent  identically-distributed  Gaussian  random 
variables  with  variance  a 2,  the  m,’s  are  the  signal  components,  N  is  the  number 
of  pixels  in  the  photodetector  array  and  /o,  A  are  respectively  the  frequency  and 
amplitude  of  the  input  signal  as  defined  in  equation  (3).  This  means  that  the  conditional 
probability  of  the  received  vector  R  given  that  the  frequency  of  the  input  signal  is  f0  is 

p(^i/4=n^^p(-(r-vu— )  (72) 


where  it  is  understood  that  the  m,’ s  are  the  corresponding  signal  components  for  an 
input  signal  of  frequency  f0  and  amplitude  A.  The  Cramer-Rao  bound  [13,  pp.  G6] 
provides  that 

Var [f0(R)~f0]  >[e 


where  f0(R )  is  any  unbiased  estimate  of  fQ. 

Taking  the  natural  logarithm  on  both  sides  of  equation  (72)  we  get 

N  ,  _  \  2 

In  (p(R  |  /„))  =  -AMn(x/2^a)  -  ^  (74) 
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where  m'  =  dm,/df0.  Squaring  both  sides  of  equation  (75)  we  get 


0\np(R  |  f0) 

df0 


-i  2 


1  N  N  N 
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.  =  1  j=  1 
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and  expanding  the  terms  of  equation  (76)  we  obtain 
2  ,  N 


d\np(R  |  /0) 


a/o 


=  ~4  ^(r.2(m!)2  +  m2(mi)2  -  2rlml(m')2 
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N  N 


+  '^(rirn',rjTn'j  —  r^m'pTijm'j  —  -f  ). 


(77) 
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Now  £(r,)  =  m,  and  E(rj)  =  rrij  and  since  r,,  r;  are  independent  provided  that 
i  ^  j,  then 

E{r,rj)  =  E(ri)E(r  j)  =  nriimj. 


Also,  it  is  easy  to  show  that  E(r2)  =  a2  +  m2.  Using  these  identities  and  simplifying 
(77)  we  get  that 

i  2 1  N  ,  i  \o 
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Combining  equations  (73)  and  (78)  we  get  that  the  Cramer- Rao  inequality  for  this 
problem  is 

Var  [f0(R)~fo\> 


E?=i  K)s 


(79) 


We  have  shown  in  section  2.0  that  for  a  sinusoidal  input  u(t )  =  Acos(2nf0t  +  d1), 

(80) 


A2 1 

mk  =  — W(/ft  -  /0) 


where  7f(/)  is  the  convolution  between  the  functions  G(f)  and  H(f) 

f°° 

«(/)  =  H(f-  f')G(f')df 

J  —oc 


(81) 


and 


II>oo 

I  w(t)exp(—i2irft)dt 
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(82) 
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As  ir  was  mentioned,  we  will  use  an  H(f)  which  is  a  rectangular  function  of  unit 
amplitude  and  width  D  Hz  symmetrical  about  /  =  0  and  a  io(t)  which  is  a  rectangular 
window  whose  duration  is  r,  the  time  taken  for  the  acoustic  wave  to  travel  across  the 
Bragg  cell.  For  this  case  we  have  that 


G(f)  = 


2sin2(27r/r/2) 

(27T/T/2)2 


(S3) 


Once  we  have  obtained  H{f)  it  becomes  easy  to  calculate  the  outputs  of  the  AOSA 
because  these  two  are  related  as  shown  in  equation  (SO).  As  an  example,  the  signal 
components  at  the  output  of  the  AOSA  are  shown  in  Figure  16  and  Figure  17.  Figure  16 
shows  the  signal  components  when  tB  =  3/4  and  Figure  17  shows  the  signal  components 
when  tB  =  6.  Both  of  these  figures  show  the  output  when  f0  corresponds  to  one  /fc  as 
well  as  when  f0  is  exactly  in  the  middle  of  two  /it’s. 

Using  numerical  differentiation  we  can  calculate  the  variance  of  the  efficient 
estimator,  Varefr  which  is  the  equality  case  of  equation  (79).  We  are  interested  in 
knowing  how  the  efficient  estimator  varies  with  tB  and  also  how  it  varies  as  a  function 
of  the  frequency  f0.  In  the  latter  case,  we  know  that  for  a  given  tB  the  efficient 
estimator  will  be  periodic  with  a  period  of  B  Hz  provided  that  the  frequency  of  the 
input  signal  does  not  correspond  to  a  frequency  near  the  edge  of  the  array. 

Obviously,  the  number  of  pixel  values  that  we  include  in  our  calculations  does  not 
have  to  be  Ar,  which  is  the  number  of  photodetectors  in  the  array.  The  reason  for  this  is 
that,  for  many  of  those  pixels,  the  values  of  the  m/s  are  relatively  small.  Hence,  we  have 
included  in  our  calculations  as  may  pixels  as  was  required  but  not  more  since  this  would 
increase  the  number  of  computations  without  any  benefit,  and  we  will  call  this  number 
n.  Using  55  pixels,  we  find  that  for  tB  =  1/4  the  efficient  estimator  is  quite  constant  as 
we  vary  the  input  frequency  and 


VVar^  =  1-27^5^ 


whereas  for  tB  =  1/2  we  find  that  it  is  also  fairly  constant  except  that  in  this  case 

\/Vareff  =  0.98 


where  x/Var^f  is  the  root  mean  squared  (RMS)  error  for  the  efficient  estimator  which  we 
will  henceforth  call  RMSeff-  Letting 
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FIGURE  16:  OUTPUT  OF  AOSA  FOR  ar  =  0,  T  =  0,  tB  =  3/4 


FREQUENCY 


FIGURE  17:  OUTPUT  OF  AOSA  FOR  ar  =  0,  T  =  0,  tB  =  6 
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FIGURE  18:  RMSeff  WHEN  rB  =  3/4 


figures  18  to  22  show  RMSeff  as  a  function  of  the  frequency  offset,  which  is  the  difference 
between  the  input  frequency  fa  and  the  corresponding  frequency  of  the  closest  lower 
pixel.  When  f0  exactly  corresponds  to  one  of  the  /*’s  then  this  frequency  offset  is  zero. 
In  these  figures  we  have  only  plotted  RMSeff  over  one  of  its  period  because,  as  it  was 
mentioned  earlier  RMSefr  as  a  function  of  frequency  is  periodic  with  period  B  Hz. 

Since  the  input  frequency  f0  could  be  any  frequency,  it  is  interesting  to  know  what 
is  the  average  RMSeff.  This  can  be  done  by  averaging  curves  such  as  those  in  Figures 
18  to  22.  This  has  been  done  for  several  values  of  the  constant  tB  and  the  result  is  as 
shown  in  Figure  23.  We  see  from  this  curve  that  the  average  RMSefr  is  relatively  small 
when  tB  is  between  0.25  and  1.25,  but  it  increases  exponentially  as  tB  is  closer  to 
2.  Figure  24  demonstrates  this  more  clearly  as  it  shows  the  average  RMSeff  for  higher 
values  of  tB.  It  should  be  noted  on  this  figure  that  the  average  RMSeff  for  values  of  the 
constant  tB  in  the  proximity  of  2,  3  and  4  are  actually  off  scale.  The  actual  value  of  the 
average  RMSeff  for  tB  =  3  is  241  /Kr  whereas  for  tB  —  2  and  4  it  is  several  orders  of 
magnitude  higher. 

Since  Figure  24  gives  us  a  lower  bound  for  any  unbiased  estimator,  it  is  clear  that 
any  such  estimator  would  have  a  very  undesirable  performance  for  rB  greater  than  1.75. 
In  fact,  we  could  even  doubt  that  an  unbiased  estimator  even  exists  for  some  of  these 
values  of  rB. 
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FIGURE  19:  RMSeff  WHEN  tB  =  1 


FIGURE  20:  RMSeff  WHEN  tB  =  1.25 
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FIGURE  21:  RMSeff  WHEN  tB  =  1.5 


FIGURE  22:  RMSeff  WHEN  tB  =  1.75 
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4.3  Peak-Detecting  Estimator 

In  tin1  previous  subsection,  we  have  presented  the  Cramer-Rao  lower  bound  on  the 
variance  of  any  unbiased  estimator  for  this  problem.  We  have  found  that  this  bound 
increases  exponentially  as  the  value  of  the  constant  tB  increases  higher  than  1.75  and 
i  hat  it  even  peaks  at  extremely  high  values  in  the  vicinity  of  certain  values  of  tD.  This 
means  that  any  unbiased  estimator  for  this  problem  would  have  a  poor  performance  for 
value.',  of  tB  greater  than  1.75  and.  in  fact  it  is  probably  impossible  to  find  an  unbiased 
i  'Umator  for  this  problem  for  some  of  these  values  of  rB. 


In  this  sul  iseetion.  in  an  effort  to  find  a  solution  to  t lit'  frequency  estimation 
-in.  we  present  the  performance  of  the  peak-detecting  estimator.  That  is.  we 
pt'-'cnt  the  performance  of  the  algorithm  that  assumes  that  the  frequency  of  the 
mpv  signal  is  the  frequency  that  •orresponds  to  the  highest  pixel  value.  Wo  want 
•<>  investigate  the  performance  of  this  estimator  since  it  is  very  natural  and  simple  to 
implement.  In  addition,  it  requires  no  modification  for  a  different  value  of  tB  or  for  a 
di::’e;vm  window  function  >c(t).  and  it  works  equally  well  on  linear  or  logarithmic  data. 


We  can  let  A’i.-Yj.  A’i . A’.v  be  A’  independent  Gaussian  random  variables  with 

I  variance  n1  for  which  the  means  are  A’j.AA.  A'3 . A'\  respectively.  This  means 

the  probability  density  functions  for  these  random  variables  are 


ft hr,  )  = 
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ami  'ince  these  random  variables  are  independent  we  can  write  their  joint  density 

tum-fion  as 

(  ,  1  IT  f  ~(x i  ~  -Y,  )2 

. ~ 
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If  we  define1  T }  as  the  probability  that  the  random  variable  X }  is  larger  than  all  the 
m.  then  it  can  be  evaluated  by  the  following  X  integrals 


Ty  =  />(  X,  <  Xj.  V/  ^  J) 


(  V  —  I  )  t i me" 


[  1  r  f  ■  ■  [  1  /(  r  i .  ,r2.  .r3 . r  v  ) 

■I  —  X  c  J  —  ~K. 

<Lr  {d.r2d.i  ;<,  ■  ■  ■  dx  }  _  i  dx  J+ j  •  ■  •  dx\  dx } 


(SC) 


htough  a  change  of  variables  can  be  rewritten 
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*,= 


'Or 


exp 


~(x  j  ~  j  )  “  \  T-r  (  1 


2n2 


n 

i=i 

>*] 


1  4-  erf 


xj  ~  X‘ 

\^2a 


dXj  ( SI 
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where 

2 

erf  (y)  =  —=  /  exp  (~t2)dt.  (88) 

\A  y0 

Using  equation  (87)  we  get  the  RMS  error  for  the  peak-detecting  estimator 

N 

RMSpk  =  *  5^(/>  -  /o)2'I'J  (89) 

N  >=1 

and  also  the  mean  error  for  this  estimator 

N 

Meanpk  =  ^(/;  -  /0)Tr  (90 ) 

i=i 

Figures  25  and  26  show  RMSpk  and  Meanpk  respectively  for  the  case  of  tB  =  0.5, 
n  =  15,  and  I\  =  A2Ir/Aa  =  10. 


FREQUENCY  OFFSET  ■  r 

FIGURE  25:  RMS  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 

( tB  =  0.5,  n  =  15,  K  =  10) 

Figures  27  and  28  show  RMSpk  and  Meanpk  respectively  for  the  same  case  except 
for  K  =  20.  Again  in  these  figures  we  have  plotted  only  one  period  of  these  functions 
since  they  are  periodic  with  period  B  Hz. 


FREQUENCY  OFFSET  r 


FIGURE  26: 


FIGURE  27: 


MEAN  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 
( tB  =  0.5,  n  =  15,  K  =  10) 


RMS  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 
(tB  =  0.5,  n  =  15,  K  =  20) 
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FIGURE  28:  MEAN  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 

( tB  =  0.5,  n  =  15,  K  =  20) 

We  see  from  figures  26  and  28  that  the  peak-detecting  estimator  is  actually  a 
biased  estimator  of  the  input  frequency.  However,  it  has  the  desired  characteristic  of 
having  a  zero  average  bias  over  many  different  input  frequencies.  This  can  be  seen  by 
the  symmetry  of  figures  26  and  28. 

As  we  have  done  with  the  Cramer-Rao  bound,  we  can  obtain  the  average  RMSp^ 
by  averaging  curves  such  as  those  of  figures  25  and  27.  This  has  been  done  for  several 
values  of  the  constant  tB  for  a  given  K.  Figure  29  shows  the  resulting  family  of 
curves  for  several  values  of  K.  These  curves  have  all  been  obtained  by  assuming 
15  photodetectors  in  our  calculations  (i.e.  n  =  15)  so  that  we  could  compare  the 
performance  from  a  common  basis.  It  should  be  noted  that,  for  any  given  value  of  K 
there  is  a  value  of  tB  which  is  really  the  smallest  that  would  be  used  in  practice.  This  is 
because  the  signal  is  so  buried  in  the  noise  that  the  corresponding  false  alarm  rate  would 
be  exceedingly  high.  We  have  not  plotted  points  beyond  that  point  on  Figure  29.  We 
see  from  that  figure  that  the  smallest  useful  tB  increases  as  we  decrease  the  value  of 
K,  which  is  in  fact  related  to  the  signal  to  noise  ratio  (i.e.  K  =  A^Ir/Aa).  We  also  see 
from  Figure  29  that  for  any  given  value  of  K,  there  is  an  optimum  tB  which  gives  the 
smallest  RMS  error.  As  in  the  case  of  the  smallest  useful  tB ,  we  see  that  this  optimum 
tB  increases  as  the  value  of  K  decreases. 

Finally,  we  note  from  Figure  29  that  for  any  value  of  K  the  average  RMS  error 
for  the  peak-detecting  estimator  asymptotically  tends  towards  a  straight  line  as  we 
increase  the  value  of  tB.  So  that  we  can  see  this  better,  we  have  plotted  a  straight  line 


VERAGE  RMS  pk  • 


tB 


FIGURE  29:  AVERAGE  RMS  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 
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on  Figure  29  that  has  a  slope  of  0.25 /t  and  that  passes  through  the  origin. 

This  behaviour  is  to  be  expected  and  to  see  why  let  us  consider  the  case  of  no  noise 
or  infinite  signal  to  noise  ratio.  Figures  30  and  31  show  the  RMS  error  and  the  mean 
error  (or  the  bias)  of  the  peak-detecting  estimator  for  this  latter  case.  It  is  easy  to  see 
that  the  average  RMS  error  in  that  case  is  B/ 4  and  hence  if  tB  =  a,  then  the  average 
RMS  error  will  be  S/4  =  a/4r  which  is  in  agreement  with  Figure  29. 


FIGURE  30:  RMS  ERROR  OF  THE  PEAK-DETECTING  ESTIMATOR 

WHEN  I<  =  oo 
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FIGURE  31:  MEAN  ERROR  (BIAS)  OF  THE  PEAK-DETECTING  ESTIMATOR  WHEN 

WHEN  K  =  oo 


5.0  ESTIMATION  OF  THE  POWER 


5.1  Introduction 

Another  important  signal  parameter  of  interest  when  monitoring  the 
electromagnetic  environment  is  the  received  signal  power.  This  information  can 
be  used  to  analyze  the  scan  pattern  of  the  transmitter  or  to  obtain  a  gross  estimate 
of  the  distance  to  the  emitter.  In  certain  system  configurations,  the  power  estimate  is 
used  for  direction  finding  and  in  any  case  we  can  know  if  the  emitter  is  approaching  or 
receding  by  comparing  its  relative  power  in  consecutive  time  frames. 

Since  we  have  no  knowledge  of  the  a  priori  probabilities  for  the  power  of  the  input 
signals,  we  are  dealing  with  a  nonrandom  parameter  estimation  problem.  This  means 
that  we  will  have  to  work  with  the  a  posteriori  density  functions. 
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In  this  section,  vve  will  also  assume  that  the  power  estimation  is  a  secondary  task  to 
the  frequency  estimation.  That  is,  we  will  assume  that  the  frequency  estimation  of  the 
signal  has  been  done  prior  to  the  estimation  of  its  power.  This  added  information  will 
allow  us  to  get  a  better  estimation  of  the  power  and  it  will  insure  that  we  associate  one 
power  estimation  with  one  frequency  estimation.  The  latter  statement  will  be  especially 
important  if,  because  of  the  system  configuration  used  the  energy  of  the  input  signal 
spreads  over  several  pixels. 

Even  though  there  are  no  general  results  on  how  to  achieve  an  optimal  estimator 
for  the  nonrandom  estimation  situation,  one  estimator  which  is  often  used  for  these  cases 
is  the  maximum  likelihood  (ML)  estimator.  This  estimator  is  based  on  the  maximum 
likelihood  principle  and  more  specifically  it  consists  in  choosing  as  our  estimate  the  value 
of  the  parameter  that  most  likely  caused  the  given  value  of  the  received  signal. 

An  important  property  of  the  maximum  likelihood  estimator  is  that  if  it  is 
unbiased  (i.e.  the  expected  value  of  the  estimate  is  equal  to  the  true  value  of  the 
parameter  of  interest)  then  it  will  satisfy  the  Cramer-Rao  lower  bound  with  an  equality 
[14].  Any  unbiased  estimate  that  satisfies  this  bound  with  an  equality  is  called  an 
efficient  estimate  because  it  is  impossible  to  find  any  other  unbiased  estimate  with  a 
lower  variance. 

In  the  next  subsection  we  present  the  ML  estimate  of  the  power  of  the  input  signal 
assuming  that  the  input  signal  frequency  f0  is  known.  We  also  analyze  the  expectation 
and  the  variance  of  this  estimate.  In  the  third  subsection  of  this  section,  we  analyze 
the  effects  of  inaccuracies  in  our  estimate  of  f„  on  this  ML  estimate  since  in  practice  f0 
would  not  be  known,  but  would  have  to  be  estimated. 


5.2  Maximum  Likelihood  Estimate 

Using  the  model  presented  in  section  2.0,  we  will  in  this  section  present  the  ML 
estimate  of  the  input  signal  power.  To  this  end,  we  let 

R  =  {r\,r2,r3,...,rN} 


be  the  received  vector  where  the  r,’s  are  the  pixel  output  values  for  a  given  frame.  We 
can  write 

r,  =  m,(/0,  A)  +  i  =  1,  2, 3, . . . ,  N  (91) 

where  the  n,’s  are  zero- mean,  independent  identically-distributed  Gaussian  random 
variables  with  variance  cr2,  the  m,’s  are  the  signal  components,  N  is  the  number  of 
pixels  in  the  photodetector  array  and  /0,  A  are  respectively  the  frequency  and  amplitude 
of  the  input  signal  as  defined  in  equation  (3). 
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We  will  assume  that  the  frequency  f0  of  the  input  signal  is  known  and  will  let 
5  =  A2  be  the  parameter  to  estimate.  This  means  that  the  conditional  probability  of  the 
received  vector  R  given  that  the  squared  amplitude  of  the  input  signal  is  S  will  be 


N  1 

p(r  i  io,  s) = n  k — y  exp 

i=l 


-(r,  -  m,)2 


2a2 


(92) 


where  it  is  understood  that  the  m,’s  are  the  corresponding  signal  components  for  an 
input  signal  of  frequency  fQ  and  amplitude  A. 


Now  to  find  the  ML  estimate  of  A2  we  must  solve 

<91n{p(£lS,/0)} 

ds 


(93) 


Substituting  equation  (92)  in  equation  (93)  and  simplifying  gives  us 

^E,=i(r.  -  mt)2/2a 2  _ 


ds 


=  o. 


Substituting  for  the  ro,’s  as  in  equation  (22)  and  simplifying  we  get 

dE!Ii(¥tt,-r,)2  _ 


dS 


=  0. 


Differ-nLatir.-  and  simnlifving  wc  obtain 


N 

E 

1=1 


(Pfs-r,ni)  =  o 


from  which  we  get  the  desired  ML  estimate 


c  4  ZL I jHj 

ML  1  ,  «? 


If  we  seek  to  find  the  mean  of  the  ML  estimate,  we  have  that 


£{^ml}  =  E 


[4  S,=i 
1 


(94) 


(95) 


(96) 


(97) 


(98) 
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which  becomes 


£{5Ml} = 


4 


N 

TntE{rt} 

i=  1 


(99) 


which  can  easily  be  calculated  to  be 

E{SMl }  =  A2  =  S.  (100) 


Hence  we  see  that  5ml  is  an  unbiased  estimate  of  the  parameter  A2.  This  is  a  very 
important  result  because  it  means  that  5ml  is  an  efficient  estimate  of  A2  and  thus  it  is 
impossible  to  find  another  unbiased  estimate  of  A2  which  has  a  smaller  variance  than 
the  ML  estimate  of  A2. 

If  we  seek  to  find  the  variance  of  5ml,  from  equation  (97)  we  have 

Vat(5ML)  =  Var^^fyr^j  (101) 


but  since  the  r,’s  are  independent  random  variables  we  have  that 

N 


Var(5ML)  = 


16 


P(Z!L ,  w?)2  U 


£«'Var(r,) 


(102) 


and  since  Var(rj)  =  <r2  for  all  i  we  get 


Var('ML)  =  nSZWY  (103) 

In  Appendix  I  we  have  tabulated  values  of  f°r  different  tB.  From 

Table  1-1  and  from  equation  (103)  we  can  see  how  the  system  design  parameters 
affect  the  variance  of  the  ML  estimate  of  the  input  signal  power.  We  see  that  the 
integration  time  of  the  photodetectors  I  and  the  Bragg  cell  aperture  time  r  are  two 
strong  parameters.  Increasing  either  of  these  parameters  will  reduce  the  variance  in  a 
square  law  relationship.  Increasing  the  value  of  tB  will  reduce  the  variance  according  to 
the  relationship  shown  in  Appendix  I.  Of  course,  if  we  reduce  the  noise  variance  we  will 
also  reduce  the  variance  of  the  ML  estimate  and  tnis  is  a  direct  linear  relationship. 


5.3  Effect  of  Inaccuracies  in  the  Frequency  Estimate 


In  the  previous  subsection  we  obtained  the  ML  estimate  of  the  input  signal  power. 
We  saw  that  the  ML  estimate  is  unbiased  in  this  case  and  this  therefore  implies  that  it  is 
an  efficient  estimate.  This  means  that  it  is  impossible  to  find  another  unbiased  estimate 
of  A2  which  has  a  smaller  variance  than  the  ML  estimate  of  A2. 

However,  these  results  were  obtained  under  the  assumption  that  the  frequency  of 
the  input  signal  fQ  was  known.  In  practice,  f0  would  not  be  known  but  would  have  to 
be  estimated.  This  means  that  there  might  be  inaccuracies  in  our  estimate  of  f0.  In  this 
subsection,  we  analyze  the  effects  of  these  inaccuracies  on  the  ML  estimate  of  the  power 
that  we  obtained  in  the  previous  subsection. 


To  this  end,  let 


s 


4  E£i  rtH, 

1  £,=i 


(104) 


be  our  estimate  of  A2  where  the  Hi's  correspond  to  frequency  fQ  which  could  be 
different  from  fQ.  Equation  (104)  is  the  same  as  the  ML  estimate  algorithm  of  equation 
(97)  except  that  it  uses  the  H,'s  instead  of  the  Hi's  to  take  into  account  the  fact  that 
our  estimate  of  the  input  signal  frequency  f0  could  be  different  from  the  actual  input 
signal  frequency  f0. 

We  are  interested  in  knowing  how  inaccuracies  in  the  input  signal  frequency  will 
affect  the  mean  of  the  estimator  defined  by  equation  (104).  To  obtain  this  information, 
we  can  take  the  mean  on  both  sides  of  this  equation 

£,l)  =  *{?&?}  (105) 


which  becomes 


N 


£<5>  =  t 

£l=l  W  1=1 


(106) 


which  after  substitution  and  simplification  becomes 


£,=i  HiHi 


£{5}  =  A‘ 


£,=i  n2 


(10?; 


Hence,  we  see  from  equation  (107)  that  the  estimator  of  the  power  defined  by 
equation  (104)  is  no  longer  unbiased  when  fQ  is  different  from  /„.  If  we  let  fa  correspond 
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to  an  jk,  which  would  be  the  case  if  we  were  using  the  peak-detecting  estimator,  then 
the  resulting  bias  would  be  as  shown  in  Figure  32,  where  we  have  assumed  rectangular 
windowing.  This  figure  shows  the  resulting  normalized  bias  for  different  values  of  tD 
which  represent  different  system  configurations.  We  have  not  included  the  curves  for 
rB  =  3/4  and  tB  —  1/2  on  that  figure  since  they  represent  very  little  variation  from  the 
curve  for  tB  =  1. 

An  important  observation  to  make  concerning  Figure  32  is  that  inaccuracies  in  the 
estimation  of  the  input  frequency  will,  on  the  average,  cause  us  to  under  estimate  the 
power  of  the  input  signal  if  we  use  the  estimator  of  equation  (104).  Small  inaccuracies  in 
the  estimation  of  the  frequency  will  have  negligible  effects,  but  inaccuracies  greater  than 
a  certain  value  will  have  considerable  effects. 

It  is  clear  from  Figure  32  that  larger  values  of  tB  present  more  robustness  to 
small  inaccuracies  in  the  estimation  of  the  frequency.  However,  as  it  was  shown  in 
section  4.0.  larger  values  of  tB  will  also  on  the  average  cause  us  to  make  larger  errors 
in  the  estimation  of  the  frequency  thereby  diminishing  the  advantage  of  larger  values 
of  tB.  In  Table  1,  we  have  tabulated  the  normalized  bias  of  the  ML  estimator  when 
the  frequency  inaccuracy  is  half  a  pixel,  for  different  values  of  tB.  We  can  see  from 
this  table  that  for  the  larger  values  of  tB  (2,3,4),  inaccuracies  in  the  estimator  of  the 
frequency  of  only  half  a  pixel  will  cause  us  on  the  average  to  under  estimate  the  power  of 
the  input  signal  by  almost  50%  if  we  use  the  estimator  of  equation  (104). 


tB 

{fo  ~  fo)r 

E{s}/A 2 

0.5 

±0.25 

0.90 

0.75 

±0.375 

0.77 

1 

±0.5 

0.63 

2 

±1 

0.54 

3 

±1.5 

0.53 

4 

±2 

0.52 

TABLE  1.  NORMALIZED  BIAS  OF  ML  ESTIMATOR  WHEN  (/0  -  fo)  =  ±B/: 


We  are  also  interested  in  knowing  how  inaccuracies  in  the  input  signal  frequency 
will  affect  the  variance  of  the  estimator  defined  by  equation  (104).  To  obtain  this 
information,  we  can  take  the  variance  on  both  sides  of  this  equation 


Var(  S )  =  Var 


(zi, «?)' 


E(s)/A2 


and  hence 


(108) 


Referring  to  Appendix  I.  we  see  that  inaccuracies  in  our  estimate  of  the  input 
signal  frequency  for  the  application  of  the  estimator  defined  by  equation  (104)  could 
either  increase  or  decrease  the  variance  of  this  estimator  depending  on  what  f0  and  f0 
are  relative  to  an  /<...  For  example,  if  we  use  the  peak-detecting  estimator  to  estimate 
the  input  signal  frequency,  then  the  variance  of  the  estimator  of  equation  (104)  can 
only  be  smaller  or  equal  to  the  variance  of  the  ML  estimator  defined  in  equation  (97). 

In  fact,  for  that  case  if  the  input  signal  frequency  is  exactly  between  two  /*’ s  then  the 
estimator  of  equation  (104)  will  have  a  smaller  variance  than  the  ML  estimator  and  this 
effect  will  be  more  pronounced  for  larger  values  of  tB. 

Hence,  we  can  conclude  that  the  estimator  of  the  input  signal  power  defined  by 
equation  (104)  performs  well  if  the  inaccuracy  in  the  estimation  of  f0  is  small.  However, 
if  this  inaccuracy  is  larger  than  a  certain  amount,  then  this  estimator  could  under 
estimate  A2  considerably.  The  variance  of  ( 104)  could  be  larger  or  smaller  than  the 
variance  of  the  ML  estimator,  but  this  effect  is  not  as  important  as  the  bias  error. 


6.0  CONCLUSIONS  AND  COMMENTS 


6.1  Summary 

We  have  presented  a  statistical  model  of  the  output  signals  from  an  acousto-optic 
spectrum  analyzer  (AOSA)  for  the  case  where  we  are  receiving  continuous- wave  (CW) 
signals.  Using  this  model,  we  have  shown  that  the  optimum  algorithm  for  the  detection 
of  a  signal  of  known  frequency  is  a  discrete  matched  filter.  In  an  effort  to  extend  this 
algorithm  to  the  case  of  the  detection  of  an  unknown  frequency,  we  have  proposed  a 
sliding  matched  filter  structure  which  can  be  implemented  as  a  finite  impulse  response 
(FIR)  filter.  We  have  investigated  the  performance  of  such  a  scheme  and  have  found  it 
to  be  quite  acceptable. 

We  also  considered  the  frequency  estimation  problem  by  calculating  the  Cramer- 
Rao  bound  which  is  a  lower  bound  for  the  variance  of  any  unbiased  estimator  of  the 
frequency.  We  saw  that  this  bound  had  very  undesirable  characteristics  for  tB  >  1.75 
therefore  indicating  the  undesirability  of  using  an  unbiased  estimator  for  these  cases.  In 
an  effort  to  find  a  solution  to  this  problem,  we  presented  the  peak-detecting  estimator 
and  we  investigated  its  performance.  We  saw  that  the  peak-detecting  estimator  is  a 
biased  estimator  although  it  has  the  desired  characteristic  of  having  a  zero  average  bias, 
and  we  saw  that  its  variance  is  relatively  well  behaved. 

Finally,  we  considered  the  problem  of  estimating  the  power  of  the  input  signal. 

We  found  it  convenient  to  approach  this  problem  by  letting  the  power  estimation  be  a 
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secondary  task  to  the  frequency  estimation.  That  is,  we  assumed  that  the  frequency 
estimation  of  a  signal  would  always  be  done  prior  to  the  power  estimation  of  that  signal. 
Under  this  assumption,  we  obtained  the  maximum  likelihood  (ML)  estimator  of  the 
input  signal  power  and  showed  that  it  was  in  fact  an  efficient  estimator,  which  means 
that  one  cannot  find  another  unbiased  estimator  with  a  lower  variance.  We  investigated 
how  errors  in  our  estimate  of  the  frequency  would  affect  the  ML  algorithm,  and  we 
found  that  it  would  always  cause  us  to  underestimate  the  power  of  the  input  signal.  For 
large  errors,  we  saw  that  the  underestimation  could  be  quite  serious.  These  errors  in  the 
estimation  of  the  frequency  were  also  found  to  affect  the  variance  of  the  ML  estimator 
but  for  the  case  when  we  are  using  the  peak-detecting  estimator,  we  saw  that  these 
errors  would  actually  reduce  the  variance  of  the  ML  algorithm. 


G. 2  Suggestions  for  Further  Research  and  Comments 

The  statistical  analyses  performed  in  this  report  were  based  on  a  system  model 
which  was  defined  in  section  2.0.  However,  since  the  field  of  acousto-optics  is  relatively 
recent,  we  are  still  in  the  process  of  discovering  new  architectures  which  can  have  certain 
advantages  depending  on  the  application.  The  system  model  for  these  new  architectures 
could  be  different  from  the  one  we  have  used  and  it  would  certainly  be  worthwhile  to 
redo  our  statistical  analyses  for  these  other  systems. 

As  we  stated  in  section  2.0,  all  our  numerical  calculations  were  done  with  the 
assumption  that  the  window  function  u’(t)  was  a  rectangular  function,  although  we 
mentioned  how  this  function  could  be  modified  to  account  for  some  other  practical 
issues  of  the  AOSA  such  as  the  attenuation  of  the  acoustic  wave  or  the  profile  of  the 
laser  beam.  It  would  certainly  be  worthwhile  to  redo  our  calculations  for  other  system 
configurations  especially  if  we  had  some  reasons  to  be  interested  in  a  specific  system 
design. 


Finally,  we  have  only  analyzed  the  case  wrhen  we  are  receiving  a  CW  signal 
although  we  mentioned  that  an  EW  receiver  would  have  to  deal  with  many  different 
types  of  signals.  Perhaps  the  most  important  type  of  signal  which  could  be  analyzed 
is  pulse  modulated  signals  including  those  with  complex  modulation  within  the  pulse. 
These  signals  are  more  complex  since  they  present  a  larger  number  of  parameters  to  be 
analyzed  each  one  of  which  possibly  exhibiting  a  large  variation  from  one  transmitter  to 
the  next.  It  is  hoped  that  this  report  has  contributed  to  a  better  understanding  of  the 
critical  issues  and  tradeoffs  that  are  involved  and  can  therefore  provide  a  basis  for  the 
development  of  an  optimal  or  suboptimal  solution  to  this  more  complex  problem. 
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In  this  Appendix  we  have  tabulated  a  few  values  of  the  expression  Tif  for 

different  values  of  tB  since  this  expression  appears  in  the  performance  calculations  for 
the  detection  in  section  3.0  and  in  the  variance  calculations  of  the  maximum  likelihood 
estimator  for  the  power  estimation  in  section  5.0.  It  should  be  noted  that  the  result 
<  of  this  summation  is  a  periodic  function  of  period  B  as  we  vary  the  input  signal 

frequency.  We  have  used  25  terms  to  obtain  these  results  since  the  remaining  terms 
were  insignificantly  small  and  we  have  assumed  rectangular  windowing.  Two  cases  have 
been  considered,  one  where  the  input  signal  frequency  /0  corresponds  to  one  of  the  jfi.  's 
of  equation  (22)  and  the  other  where  /0  is  exactly  between  two  fk s.  The  former  case 
represents  a  maximum  for  a  given  value  of  rB  while  the  latter  represents  a  minimum. 
The  results  of  these  calculations  are  shown  in  Table  1-1. 


tB 

Ei, 

fo  corresponds 
to  one  fk 

fo  between 
two  /*’  s 

0.25 

0.163 

0.163 

0.50 

0.308 

0.308 

0.75 

0.452 

0.396 

1.00 

0.611 

0.409 

1.25 

0.730 

0.415 

1.50 

0.792 

0.434 

1.75 

0.814 

0.449 

2.00 

0.817 

0.452 

3.00 

0.868 

0.467 

4.00 

0.903 

0.475 

5.00 

0.920 

0.480 

6.00 

0.934 

0.483 

TABLE  1-1: 
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